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Abstract 



The combined capacity expansion, production- 
inventory problem Is modeled as a linear. Integer 
program and then extended to also allow capacity de- 
crease. The models assume constant returns to scale 
In the production function of a firm which must meet 
deterministic demands for a single product over N 
periods, with no backordering, at minimum cost. 

A linear transformation Is used to obtain 
equivalent forms of the models which are then decompos- 
ed Into fixed cost and variable cost parts. Enumer- 
ation Is done on the fixed cost variables and transpor- 
tation sub-problems are solved with the remaining var- 
iables. Special demand and cost structures are Inves- 
tigated, one of which leads to formulation of the en- 
tire problem as a network. 

Extensions are made to Include plecewlse- 
llnear, convex production costs and lumpiness In ca- 
pacity acquisition. A reinterpretation of the models 
to handle the employment smoothing problem Is also pre- 
sented. Finally, computational results of a Fortran IV 
coding of the basic algorithms are presented. 
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Chapter I 



I.l Introduction 



The quantitative approach to the solution 
of problems in the realm of managerial economics can 
be traced back to the efforts of the early Industrial 
engineers in the United States at the end of the nine- 
teenth century. These men, most notably Frederick W. 
Taylor, pioneered the Industrial engineering concepts 
such as Job description, time-motion studies, etc., 
which had such a great influence on the development 
of the assembly line in Industry during the first 
half of the twentieth century. 

With World War II, however, an entirely new 
set of tools was refined and applied by the early prac- 
titioners in the Infant field of Operations Research. 
Though most of the techniques had their origin in 
military applications, they gradually found their way 
into the industrial sphere in the peace-time economic 
expansion resulting from the war. In the intervening 
decades , the complexity of the practical problems 
solved has moved apace with the rapid advances in theory 
typical of a new field. The logical result is that am 
increasing amount of effort has been put into Inte- 
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gration of past techniques using newer theory in an 
attempt to provide more realistic approaches to problems 
with broader horizons. 

The present thesis makes such an attempt by 
offering models of the combined production-inventory 
and capacity expansion problem and algorithms which 
give all Integer solutions to these models. Dis- 
cussion of the importance of the combined problem and 
the Integral approach are given later in this chapter. 
This is the first Integer approach to the combined 
problem known to exist. 

The basic technique used will be branch and 
bound with network sub-problems. A unique feature 
of the algorithm is the capability of starting it 
at any point in the enumeration. This is especially 
meaningful when linear programming is used to give a 
universal, and often non-integer, optimal solution. 

The enumeration can be started from a suitable point 
with respect to the linear programming optimum to 
determine the existence of any optimal, rounded-off 
solutions rather quickly. 

The famous economic problem of ” lumpiness” 
in capacity acquisition as well as multi-shift, and 
overtime, operation of the production facility will 
be treated. A discussion of the relationship between 
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the traditional capacity expansion and employment 
smoothing literatures and re-interpretation of the re- 
sults in terms of employment smoothing will also be 
given. 

1.2 Problem Description 

A useful paradigm is a company with a single 
producing facility and a homogeneous product faced 
with a schedule of deterministic demands over a set of 
discrete periods to some pre-set horizon. These de- 
mands must be met on time from production or existing 
inventory at minimum possible cost. 

The capacity of the production facility is 
defined as its output in units per time period when 
equipment is being run at its maximum efficient rate 
with a single shift not working overtime. The models 
will be extended to allow multiple shifts and over- 
time work, but note that these additions are not ger- 
mane to the basic definition of capacity. 

The company's basic problem, assuming that 
cumulative demand never exceeds cumulative capacity 
(where cumulative capacity now includes all regular 
and overtime production), is to find a production- 
inventory schedule to meet demand at minimum 
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cost,^ References to this problem are ublquitlous in 
both the operations research literature and the ec- 
onomics literature, and will not be discussed in this 
work. 

If, however, cumulative demand does exceed 
cumulative capacity within the planning horizon, then 
the company must make a capacity expansion decision. 

That is, it must decide at what time Cor times) and by 
what amount Cor amounts) to invest in capital equipment 
so that production can satisfy demand. The intuitive 
difficulty of this problem is the fact that a change 
in the capacity in one period affects the production- 
inventory schedule for all future periods, thus lending 
combinatorial aspects to the capacity expansion problem. 
Further difficulty arises when capacity must be added 
in indivisible units, each capable of producing large 
numbers of units of output as is usually the case in 
industry. 

Capital equipment itself has many aspects. 
Deterioration might be a significant factor with the 



^The existence of capacity in excess of demand 
does not preclude a decision to expand capacity since 
there could exist circumstances where such a move 
would lead to lower costs. 
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type of capital Installed, in which case the capacity 
expansion problem Is further complicated by the need 
for equipment replacement. On the other hand, fixed 
capital might have a relatively long life, making per- 
Iodic maintenance the only slgnlflcajit factor within 
the time-span of Interest, It Is also possible that the 
fixed capital equipment can be salvaged or scrapped If 
warranted by changes In the demand situation. The first 
case can only be approximated by the models which 
follow and then only by allowing maintenance costs to 
reflect the Increased Investment over time needed to 
maintain equipment at Its original capacity. In add- 
ition to not obtaining equipment replacement decisions 
from the models, the a^>proxlmat Ion Is faulty In that 
the Increased maintenance costs are applied uniformly 
to equipment In each period, regardless of age. The 
latter two cases are modeled exactly, however, the sec- 
ond by the explicit structure of the maintenance costs 
and the last Implicitly by the relevant mathematical 
programming model. 

I . 3 Economic Background 



With respect to the economic geography, the 
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problem lies on the boundary of capital theory and the 
classical theory of production. A good discussion of 
the problems arising from the merger of these two dis- 
ciplines is contained in the unpublished Ph. D. thesis 
of Carroll (Carroll 68). 

Apropos of capital theory, a decision must 
be made as to the proper points in time, and the proper 
amounts, for Investment in capacity; however, there is 
no requirement present for alternative Investments. 

In general, if expansion of capacity in a number of 
different areas were being considered in each period 
with a limited budget, then the extensive work in cap- 
ital budgeting would be relevent.^ Further, capital 
budgeting problems typically abstract from the micro 
problem of optimizing the form of each investment pro- 
ject, being concerned mainly with the macro allocation 
of funds to projects. The micro-optimization is an im- 
portant part of the work that follows. Extensions of 
the present problem to the capital budgeting sphere 
will be discussed in the final chapter. 



^A classical application of mathematical program- 
ming to the capital budgeting problem is given by 
Weingartner (Welngartner 63 ). 
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In the theory of production as dlelcussed by 
Carroll, consideration of production, inventory, and 
capacity expansion places the problem in between the 
short term, where capacity is fixed and decisions are 
only made on levels of production, and the long-run de- 
cision of optimal overall capacity level, which ab- 
stracts from all short-run decisions. In effect, the 
mechanism of inventorying production functions as the 
mediator between the short-run minimum cost production 
problems and the long-run minimum cost capacity expan- 
sion problems, serving to "smooth” the extremes in the 
optimal solutions to the individual problems. 

A more satisfying concept of the spectrum of 
decision time-spans is given by Maxwell (Maxwell 65). 

He considers the very short-run as the time-span in 
which all factors of production are fixed in both ac- 
quisition and use. Successive decision time-spans are 
encountered as factors become variable in use and ac- 
quisition, giving a succession of short runs of in- 
creasing length. The long-run is finally achieved when 
all factors become variable in both use and acquisition. 
Against this spectrum, the problem at hcind ranges from 
a succession of relatively long short-runs when only 
capital stock is fixed in acquisition and all factors 
are variable in use, to a long— inan over which the 
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capital stock is variable in acquisition at a succession 
of discrete points. Correspondence with classical 
theory is not exact though, since capital is durable 
and not perfectly flexible as required in the classical 
long-run. The necessity for providing inter-period rela- 
tionships through the mechanism of Inventorying in the 
case of Inflexible capital was pointed out by Carroll: 



...So long as the theory of production ab- 
stracted from the particular problems of fixed 
in-puts, it was essentially a mono-period 
theory in which inter-period relationships 
could largely be neglected. Even the long- 
run analyses tended to treat the question of 
determining the optimal capital stock as a 
once-and-for-all decision. Since capital was 
generally assumed to be completely flexible 
in the long-run, questions of the impact of 
one period's capital accumulation decision 
upon that of another were not raised. But 
once capital came to be viewed as a durable 
good, which was not perfectly flexible, 
inter-period (l.e., between long-run periods) 
relationships were a matter of Interest. 

(1968, p. 5) 



In the following work, a salvage value of capital 
goods is introduced and constrained to be less than the 
acquisition cost, thus ensuring at least a mild flex- 
ibility in capital. 

Maxwell also presents a variation on the con- 
ventional theory of the firm which is relevant here. 

The contrast is in the derivation of the production 
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function which, in the classical theory, is obtained 
from the combination of varying Inputs of consumptive 
factors of production with the fixed factor to obtain 
a graph of maximum outputs over the range of the con- 
sumptive factors. The underlying assumption is, of 
course, that the fixed factor is indivisible in use and 
must be used to the available limit. If Instead, the 
fixed factor is viewed as indivisible in acquisition, 
but divisible in use, a more useful approach to the 
firm's production function can be derived. 

Assume constant returns to scale in the pro- 
duction function. Then, for each level of input of 
consumptive factors, the fixed factor is used up to its 
zero level of marginal productivity giving a linear ex- 
pansion path for the firm up to the capacity of the 
fixed factor. Thus, under constant returns, the scale 
of use of the fixed and consumptive factors of produc- 
tion is varied linearly up to the capacity of the fix- 
ed factor, while the proportions of both remain cons- 
tant. Beyond the capacity of the fixed factor the pro- 
portion of the consumptive factor is increased, giving 
returns to proportion in production as opposed to the 
returns to scale obtained in the former range. The 
importance of this concept to the theory of the firm is 
pointed out by Maxwell; 
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Many empirical cost studies Indicate that the 
marginal costs tend to be constant very nearly 
up to capacity, suggesting that this range In 
which less than all of the stock of fixed- 
factor services Is being utilized may Indeed 
be a slgnlflcamt portion of the total range 
of outputs that the firm Is capable of pro- 
ducing In the short-run. (1965, p. lo6) 

For the problem of Interest, the range of the 
production function reflecting returns to proportion Is 
of no Interest since one of the decision variables for 
each period Is the optimal level of fixed capital for 
that period. Thus, concern can be focused solely on the 
constant returns to scale portion of the production 
function assumed present up to the capacity of the fix- 
ed factor. 

In the problem which follows, the bundle of 
consumptive factors of production 9^ which Is needed 
In period 1 for one unit of output Is assumed known. 
Further, It Is assumed that the fixed factor of pro- 
duction Is used at Its zero marginal productivity level 
for this bundle 9^^ or any multiple of It. The capacity 
of the fixed factor In any period can then also be 
expressed In terms of the maximum number of bundles 9^^ 
which can be used efficiently with the fixed factor In 
the period (In this case the period Includes only reg- 
ular time work by one shift). The level of production 
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in period i which Is must thus satisfy 

1 = 1, . . . ,N 

These assumptions, when taken In the context 
of fixed capital Investment, suppress the problem of 
labor force. That Is, although the labor force Is a 
non-consumptive factor of production, the necessary 
assumption here Is that labor force la always adjusted 
to the exact size required to run the capital equipment 
at the predetermined optimal level Xj^ In period 1. An 
Interpretation of the problem of capacity In terms of 
labor force rather than capital equipment will be 
given subsequently. 

The objective Is to minimize the present 
value of the costs Incurred over N periods to meet a 
predetermined demand schedule. It should be noticed 
that the problem Is abstracted from pricing consider- 
ations and from determination of optimum output 
schedules; In fact, concern Is with the optimum method 
of meeting any demand schedule. To achieve the pre- 
sent value criterion. It Is only necessary that the 
Individual costs In the objective function of the model 
be discounted by the appropriate cost of capital In 
period 1. Thus discounting does not have to be modeled 
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but can be introduced explicitly in the relevant model 
coefficients. The subject of the proper choice of the 
discount factor p is not germane to this thesis and 
will not be discussed. 

The major stumbling block to the merger of 
the capital and production theories has historically 
been the problem of fixed capital as a non-consumptive 
factor of production. This question is also lucidly 
discussed by Carroll. To summarize, some method must 
be found to measure the value of a factor of production, 
i.e* fixed capital, which is not consumed in the pro- 
ductive process, but which must be present for produc- 
tion to take place. In this thesis, the value measure 
will be the reduction in cost obtained in the produc- 
tion of a given demand schedule, a measure which is im- 
plicit in mathematical programming optimization. 

Finally in the economic context, mention 
should be made of the historic bugbear of marginal 
analysis. Indivisibility in acquisition of capital 
stock. This discreteness in acquisition possibilities 
destroys the differentiability of most useful functions 
at best, and their continuity itself at worst. As will 
be described later, a simple adjustment in the integer 
programming formulation will hcindle this problem 
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fairly easily and, In some cases, even improve the al- 
gorithms. 

Problem Formulation 



The Interval of time under consideration is 
divided into n periods of equal lengths with the first 
period designated as 1 and the n^h^ ©r horizon, period 
designated N. It is assumed that demand in period i, 
dj^, is known deterministically for i = 1 ,...,n. Then 
using the following variables: 

production in period 1, 1 = 1,...,N 
v^- inventory from period 1 to period 1+1 
Wj^- capacity added in period 1, 1 = 1,...,N 
the following coefficients; 

Pj^- cost of a unit bundle S^^of consumptive 
factors of production in period i. 



hi- 



c<- 



mi- 



1 « 1 , . . . ,N 

unit Inventory cost from period 1 to 1 + 1 
construction cost of a unit of capacity 
in period i 

maintenance cost of a unit of capacity in 
period 1 

where coefficients are assumed integer and 
p^ — 0, h^ ~ — 0, mj^ — 0 1=*1,...,N 



> 
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and the following explicit Integer variables: 

Wq- Initial capacity (w^ - 0) 

Vq- Initial Inventory (v^ - 0) 

If all functions are assumed linear, the capacity ex- 
panslon-productlon-lnventory problem (CPI) can be 
written as the following Integer programming problem; 



PI: 



min 



N 1 

^2^CPlXi+hiVi+ciWi+mi* wj+w^p] 



s. t. 

(Pl.l) 



^1 



1 

I w. 
J = 1 J 




1 - 1 , 



,N 



Xl - Vl - dj^ - Vq 

(PI. 2) x^ + - Vi - di 1 - 2,...,N-1 

""N + '^N-l “ 

X, V, w - 0 and INTEGER 



where the right-hand side of the constraints Is assumed 
Integer. Constraints (Pl.l) restrict production to be 
less than or equal to capacity In any period, and con- 
straints (PI. 2) require demand to be met exactly out 
of the production and Inventory In each period. All 
costs are non-statlonary and are explicitly discounted 
to the first period by the appropriate discount rate p^, 
which Is external to the model as discussed In section 
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1.3. 

Initial inventory, Vq, and initial capacity, 
Wq, are set explicitly to the values relevant to the 
problem being solved. Initial Inventory and initial 
capacity are rather straightf onward in the problem con- 
text and no difficulty in obtaining them should occur 
in the real world. Pinal Inventory, however, is more 
of a subjective variable. If Vj^ is allowed to be an 
implicit variable, then Vj^ 0, since if v^^ > 0 then 
the level of production in periods N, N-1, ... could 
be reduced until Vj^ ■ 0 , giving a lower cost while 
still minimizing cost (recall that cost coefficients 
are non-negative). On the other hand, Vj^ can be use- 
ful as an explicit variable as shown by an example. 
Assume, for N » 3» that the demand schedule is 
d*^ « (1,7,20).^ In this case, Vj^ could subjectively 
be set to some value greater than zero as a surrogate 
for an expectation of continued Increase in demand be- 
yond the horizon. 

The objective function in PI can be re- 



3ln the remainder of this thesis, an untransposed 



vector is assumed to be a column vector 
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written as 

min Z-i “ Z [pi x<+hi vi + (c4+ ^ m.)w. ] + Wq» ^ m. 

■L i-1 ^ ^ ^ ^ J-i J ^ 1-1 i 

It can be clearly seen from this aggregated function, 
that the total cost of an incremental unit of capacity 
is equal to its acquisition cost in period i plus 
the discounted sum of the malntencince costs for one 
unit of capacity up to the horizon. Maintenance costs 
themselves are period dependent ( l.e.. Independent of 
the age of the Installed capacity). 

The CPI model can be extended easily to a 
capacity increase-decrease, production-inventory model 
(CIDPI) as shown below. Let: 

rj^— capacity scrapped in period 1 
fj^- salvage value per unit of capacity 
scrapped in period i 

all other variables and coefficients are as before and 
assume 

f^ ^ 1 - 1, ... ,N 

fj + ^ J > 1, 1 = 1,...,N-1 

Then the CIDPI can be written as the following Integer 
programming problem: 




% 

p 



P2: 



min Z 2 “ 



17 . 



N r 

^Z^CPiXi+hiV^+CiWi-f^rj^+mi* ( 



J 




s .t . 

1 

(P2.1) - L (wj-rj) - Wq 1 = 

J ™ 1 

XI - Vi » dl - Vq 

(P2.2) x^ + vi_i - ’'^1 “ 1 ■= 2,...,N-1 

XN + 

i > 

(P2.3) Wq + E (wj-rj ) - 0 1 » 1,...,N 

J “ 1 

X, V, w, ?-0 and INTEGER 



where the right-hand sides of the constraints are in- 
teger as before, and constraints (P2.1) and (P2.2) have 
the same interpretation as previously. Constraints 
(P2.3) ensure that current capacity is always non- 
negative. 

It should be noted here that P2 contains 
the flexibility of allowing f^ < 0 as might be the case 
if, instead of a salvage value for equipment, there 
were a positive cost of scrapping capacity in some 
period (in which case scrapping might be a means of 
avoiding unnecessary maintenance cost on unneeded 
equipment). For f^^ - 0, P2 is bounded because of the 
initial assumptions made on c^^ and f^, and for fj^ < 0, 
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P2 is bounded because of the form of the objective 
function. 

As with PI , the objective function of P2 
can be aggregated as follows: 



min Z2 



N N 

^^^CPlXi+hiVi+(Ci+^E^mj)wi 



+ 



w. 



N 

• Z 

1=1 



mi 






PI becomes a special case of P2 if the 
scrapping cost Is set to Infinity (i.e., fi= -QOVl) in 
P2 . Then r = 0 always and constraints (P2.3) become 
redundant, giving back the form of PI . Two distinct 
models will be retained, however, for both economic and 
theoretical reasons. 

Finally, in the CIDPI model, the value of 
fixed capital is central to the minimum cost solution. 
During periods 1 to N, this value is reflected in the 
salvage coefficient f. To be complete, though, the 
model must include the value of fixed capital after 
period N, and here two cases are relevant; either the 
firm expects to go out of business after period N or it 
expects to continue production, but its planning horizon 
is limited by dependable demand projections. In both 
cases, a dummy period N + 1 can be added to the model 
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with dj^+i = 0, Then, in the former case, fN+1 is Just 
the relevant salvage value of the fixed capital, while 
in the latter case, equal some dis- 

counted estimate of the future value of production 
possible with the Installed capacity. 

As for the CPI model, nothing need be done 
since the salvage value of capital has no bearing on 
the problem. That is, capital only has value in the 
problem context to the extent that it can meet the de- 
mands. 

I. 5 Related Literature 

A large literature exists in the area of pro- 
duction planning and there is at least a modest one on 
capacity expansion. But the combined problem, though 
economically significant, has been given relatively 
little attention, and no attention at all has been 
given to the necessity of providing good Integer solu- 
tions to such a problem. A related literature closely 
allied with the combined problem does exist, however, 
in the area of employment, or work-force, smoothing. 
Each of these areas will be discussed subsequently. 

Mention must first be made of the relation- 
ship among the existing capacity expansion literature 
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economics, and reality. Any demand-related problem 
naturally depends on a reasonable projection of demand, 
either stochastic or deterministic, within a believable 
horizon. Since forecasting is, at best, a risky busi- 
ness, the veracity of any demand projection must have 
an inverse relationship to the length of the horizon 
over which the projection is made, A reasonable figure 
for most companies would probably not exceed five years. 
Within such a time span, however, costs of production 
and inventory become economically significant, while 
beyond this span, any quantitative approach rapidly 
approaches speculation. The only valid conclusion then 
would seem to be that a realistic, quantitative 
approach to any capacity problem in the production sec- 
tor must Include the effects of the production and 
Inventory policies. 

A further economic and realistic condider- 
ation in the capacity expansion problem must be the 
fact that capacity, in whatever form it may take, is 
usually obtainable only in fixed Increments. Further- 
more, an incremental acquisition of capacity could, by 
the nature of the machine or whatever that is added, 
correspond to a large incremental capability in terms 
of units of output. Thus, although the roundoff of an 
exact, non-integer solution may give cin answer that 
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Is "close" to optimal in terms of the demand profile, 
the error in terms of cost could be very large. It 
must be concluded then that, at least in the context of 
capacity expansion, integer solutions are important. 

Perhaps the classic early delineation of the 
capacity expansion problem Has presented by Manne 
(Manne 6l) in his analysis of the basic problem for 
both deterministic and stochastic demands. Early 
linear programming models of the problem were pre- 
sented by Fetter (Fetter 6l) and Veinott and Wagner 
(Velnott-Wagner 62). An equivalent form of the 
lattera* model and a corresponding network formulation 
can easily be obtained from P2 by setting production 
and Inventory costs to zero, Kalro and Arora (Kalro- 
Arora 70) extended the basic linear programming form- 
ulation to handle stochastic demands. 

A dynamic programming approach to capacity 
expansion is given by Howard and Nemhauser (Howard- 
Nemhauser 68) for deterministic demands and extended 
to the stochastic demand case by Montgomery (Montgomery 
71). An extension of the dynamic programming model 
for deterministic demands to handle concave costs of 
expansion is made by Manne and Veinott (Manne- Veinott 
67). None of the above work takes into account the 
production-inventory problem, however, and thus will 
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not be considered as directly antecedent to this thesis. 

Zabel (Zabel 63) was one of the early authors 
to look at the total problem, but his assumption of no 
Inventory effectively reduces the problem he solves to 
one of capacity expansion. Hlnomoto (Hlnomoto 65) 
essentially retained Zabel’ s no-inventory assumption 
and extended his work to include technological Impro- 
vement, thouch he couches the work in tenns of replace- 
ment theory (falling to mention Zabel). 

Iglehart (Iglehart 65) seems to have analyzed 
the first complete model for the CPI problem. He ex- 
tends Zabel ’s work to allow inventory and restricted 
back-ordering using analytical techniques to establish 
optimal capital accumulation and production policies. 

He poses what seems to be a serious restriction, how- 
ever, when he assumes that the consumptive factors of 
production are always spread evenly over all fixed cap- 
ital, in effect eliminating the possibility of idle 
capacity. Carroll (Carroll 68) extended Iglehart' s 
work to stochastic demands and pricing considerations 
using dynamic and non-linear programming techniques to 
establish optimal capital investment and production 
policies. 

The most significant general formulation of 
the CPI problem to date, from a quantitative viewpoint, 
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has been the work of Coblan and Lele (Coblan-Lele 70). 
They give a general dynamic programming formulation of 
a multi-plant, multi-inventory facility firm faced 
with deterministic demands and show how the model can 
be hierarchically decomposed to simplify the solution 
process. Unfortunately, even when only one plant and 
one Inventory facility are assumed, the model still 
contains two state variables and gives non-integer solu- 
tions . 

The area of recent work which appears most 
relevant to the area of capacity expansion seems to 
have been done in the related literature on employment 
smoothing. Interpretationally , work force level can 
be viewed as a measure of capacity; the assumption 
made in this literature being that the fixed capacity 
level is given in each period, as opposed to assuming 
that the work force itself is optimally adjusted to 
the level of output, as mentioned in the previous sec- 
tion on economics. The greatest objection overcome in 
the employment smoothing literature is that, by nature 
of the problem, the production-inventory policy is 
always Included. To date, however, integrality has not 
been considered an important aspect of the problem. 

An excellent delineation of the problem and methods of 
solution plus a bibliography through 196? is given by 



Silver (Silver 67). 

A linear programming model which is close in 
form to P2 and a heuristic for "rounding-off” the 
linear programming solution are given by O'Malley, 
Elmaghraby, and Jeske (O'Malley et al. 66). As is well 
known, and demonstrated by example later in this thesis, 
rounding-off can be far from optimal and a very costly 
approximation to the optimal integer solution. 

The major difference between P2 and the 
model of O'Malley et al, is the non-statlonarlty of all 
costs in the former, as opposed to the complete station- 
arity of costs in the latter. In addition, the capa- 
city decrease coefficient in P2 can be positive or 
negative, within the bounds of the preconditions, to 
correspond to salvage value or scrapping cost re- 
spectively, while the corresponding term in the O'Malley 
model is a firing cost and constrained to be positive. 
Lastly, P2 is an integer program as opposed to the 
linear program of O'Malley. Still, O'Malley's model is 
not a special case of P2 since it has a term for over- 
time production. This, however, is a simple extension 
of P2 done in Chapter III, at which time the O'Malley 
model will become a special case (except for integral- 
ity). 

The balance of the recent major work in 
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employment smoothing has been done in two companion 
papers by Lippman, Rolfe, Wagner, and Yuan (Lippman et 
al. 6 ?a and 67b) and a later paper by Lippman and 
Yuan (Lippman-Yuan 69). The model addressed in the I967 
papers is a generalization of the linear model of 
O'Malley et al. with convex-like, non-stationary pro- 
duction costs, increasing and non-stationary inventory 
costs, and v-shaped, stationary smoothing costs. De- 
mand-dependent upper and lower bounds on work force 
sixe are given find optimal policies are developed for 
the special cases of monotone increasing and monotone 
decreasing demands. Algorithms are then given for the 
special cases of monotonicity and the added restriction 
that all costs are linear, though non-statlonarity is 
retained for production and Inventory costs. 

In the 1969 paper, similar analytical work 
is done on a model which is the same as that of 1967 
with the exception that smoothing cost is permitted to 
be non-stationary. This paper also points out that the 
upper bounds and algorithms developed in the I967 
papers are invalidated under the more generalized 
smoothing costs. No further algorithms are presented. 

The models in this thesis are special cases 
of the work of Lippman et al. only Insofar as the more 
general form of the latter's production and inventory 
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cost functions are concerned. Integrality and the 
existence of algorithms In this work seem to make 
further comparison unnecessary. The most cogent point 
In the above discussion would seem to be the recogni- 
tion of the Interrelationship between capacity expaui— 
slon literature and employment smoothing literature, 
about which more will be said In the final chapter. 

1,6 Plan of the Work 

In Chapter II, PI and P2 will be trans- 
formed to the equivalent problems on which the algor- 
ithms given In the chapter are based. Bounds are de- 
veloped here for the variable which acts as a surrogate 
for capacity. Finally, more efficient bounds and 
algorithms are presented for special monotone demand 
and cost cases. 

Chapter III presents a network formulation 
of the entire problem for the special case when cap- 
acity Increase eind decrease costs are equal In each 
period, but are allowed to be non-statlonary . The 
basic models are then extended to Include Indivisibil- 
ity In acquisition of capacity and plecewlse-llnear , 
convex production costs that provide for consideration 
of multi-shift work and overtime. 
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A discussion of similar Integer programming 
algorithms and a reinterpretation of the results of 
this thesis for the employment smoothing problem are 
presented In Chapter IV. Also presented here are the 
results of test problems run with a FORTRAN IV coding 
of the algorithm along with a summary and conclusions. 



I 



Chapter II 



II. 1 Introduction 



Models PI and P2 of Chapter I could be re- 
formulated to eliminate the Inventory variable v, in 
which case PI would have 2N constraints and 2N variables 
and P2 would have 3N constraints and 3N variables. In 
the capacity expansion context with a period of one 
year and N = 5, neither problem would be very large and 
existing all integer algorithms could be used. On the 
other hand. In the employment smoothing context with a 
period of one month and N = 24, where furthermore the 
CIDPI model Is most relevant, the problem size is out 
of the range of all except branch and bound type al- 
gorithms. 

Use of implicit enumeration seems out of the 
question since the variables of PI and P2 can take on 
rather large values, and expansion to 0 - 1 variables 
would quickly exceed the capacity of existing algorithms 
for even small values of N. Primal and dual cutting- 
plane algorithms have the disadvantages of highly un- 
predictable convergence and, for the latter, lack of a 
current feasible solution until optimality is reached. 
Although there is the possibility that one of these 
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other methods possesses an extremely efficient adapta- 
tion for this special problem, this thesis Investigated 
the most obvious approach- branch and bound. In this 
chapter, special branch and bound algorithms will be 
developed for PI and P2 which exploit the structure In- 
herent In these problems. 

In section II. 2, a linear transformation Is 
done on PI to obtain an equivalent form P3 for which an 
algorithm can be more easily developed. The nature of 
the algorithm is discussed subsequent to the transform- 
ation. When treated as a function of certain of its 
variables, the optimal value of the objective function 
of P3 has convexity properties over certain ranges of 
those variables which are established in section II. 3 
by a series of lemmas, theorems and corollaries. The 
properties of the integer program P3 plus the results 
from section II. 3 lead to the set of bounds and op- 
timality tests given in section II. ^ which are sub- 
sequently used in the branch and bound algorithm of 
section II. 5. 

A similar development Is followed for the 
CIDPI model P2 in sections II. 7 and II. 8. The linear 
transformation used Is the same and the resulting 
model P5 has properties similar to those of P3 both In 
problem structure and objective function convexity. 



I 
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The modifications to the basic algorithm of section 
II. 5 needed for P5 are given at the end of section II. 8, 
For both transformed models there are special 
cases of demand and cost structures which lead to major 
simplifications in the algorithms. These structures, 
their properties, and the attendent modifications of 
the algorithms are given in section II. 6 for the CPI 
model and section II. 9 for the CIDPI model. 



II. 2 Transformation of PI 



A set of necessary and sufficient conditions 
for the feasibility of a w can be obtained from con- 
straint sets (Pl.l) and (PI. 2). The method in general 
is to solve the inequalities (Pl.l) for x by adding 
the appropriate terms to each side. Then for any 
k = 1 , 2 , . . . ,N , sum the first k constraints of (PI. 2) to 
obtain a constraint of the form 

k 

+ X2 + ••• + - ''k ■ - "o 



and substitute for the from (Pl.l) giving an in- 
equality of the form 



kw^ + (k-l)w2 + 



"'k - 




kw 



o 



• • • 
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Repeating for each value of k gives the set 



Wi - Vi - di - Vo - 



( 2 . 1 ) 



2wi + W2 - V2 - + d2 - Vq - 2w^ 



+ ... + Wj^ — Vj^ — dj^ + ... + dj^ -> Vq — Nw^ 



If a vector of surplus variables t - 0 Is 
added to (2.1) and the v^^ and t^^ In each constraint of 
(2.1) are aggregated as 



Yl = Vi + tf 1 » 1,. . . ,N 

then y - 0 since v - 0 and t - 0, and (2.1) can be 
written In vector form as 



( 2 . 2 ) 

where 



Bw - jr = 




w 



o 




0 

1 

2 



0 

0 

1 



...00 

...00 

...00 



N (N-1) (N-2) ...21 



and the form of the vectors is obvious 
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Solving for w from (2,2) gives 



/s 

w ■ B"^d + B”^y - B“1 vq - B“1wq 



where 




1 

-2 

1 

0 



0 

1 

-2 

1 



0 0 
0 0 
1 0 
-2 1 



0 

0 

0 

0 



0 0 
0 0 
0 0 
0 0 



0 0 0 0 ... 1 -2 1 



Expansion gives the following set of relations 

Wi = di + yi - Vq - w^ 

W2 - d2 - di + y2 - 2y^ + v^ 

(2.3) 

W3 = d^ - d2 + y3 - 2y2 + y^ 

“n ’ - %-i * >'n - 2yN-i ♦ yN-2 

which expresses w explicitly In terms of demand. 
Initial Inventory, Initial capacity, and the newly In- 
troduced aggregate variables y. The structure of the 
right-hand side Is rather unique and, as will be seen 
shortly, extremely useful. The above transformation 
will be referred to as Standard Transformation 1. 
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An Interesting economic Interpretation can be 
given to the variable y. If overcapacity is defined as 
the cumulative excess of capacity over cumulative de- 
mand up to and Including period 1, then y^ is the over- 
capacity in period 1. For instance, ifN = 3, d'^*(3, 

5 , 6 ), and * ( 5 , 0 , 0 ), and = v^ = 0 then y]^ * y 2 ® 2 
and y^ = I 4 Further, in any period 1, v^ is the por- 
tion of overcapacity y^ which is used to produce for 
demand after period 1 and stocked as inventory, and t^^ 
is the cumulative amount of idle overcapacity as of the 
end of period 1, Thus, if* = 0, then yjq is the total 
idle capacity in the model. 

Model PI can be transformed by substituting 
for from Standard Transformation 1 to give the 
following equivalent model. ^ 



P3: 



min Z- 



i=l 



[p 



1^1 + + =;yi^ + 



K 



3 * t • 



(P3.1) 



Xi < di + yi - v^ 

Xi - di + yi - y^_j^ 1 « 2,...,N 



^For the remainder of this thesis, the objective 
function of P^ will be referred to as 



3 ^. 



Xi - - di - 



(P3.2) 









^1 - '^o + '^o - ^1 
(P3.3) 72 - 2y^ - ^1 - ^2 ~ '^o 

y± - 2yi.i + y±_2 ~ ^±.i “ ^i i = 3 ,...,n 
3c, V, y - 0 and INTEGER 



where constraint set (P3.3) comes from substitution in 
w ~ 0 in PI; the restriction of y to Integers comes from 
the restriction of w to Integers in PI and the fact 
that the right-hand side of constraint set (P3.3) is in- 
teger by assumption, and where 



c^ » Cj^ - 2cj^^^ + Cj^^_2 + m^^ - i “ l,...,N-2 



N-1 



Cn- 1 - 2cj^ + - mjj 



"n “ "N ^ \ 



N-1 






If, in P3, a feasible y is given then the problem 
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becomes 

N 

P^: min = Z 



S • b • 



(P4.1) 



(P4.2) 



- di + yi - 



- di + yi - yi_i 1 = 2,...,N 



Xt - V. = d, - V 

1 1 1 o 



H * '’1-1 - ''i =■ “i 



1 = 2 N-1 



‘n ''N-1 ' * ''n 



X, V - 0 and INTEGER 



Which Just Is the production-inventory problem, para- 
metric in y; but which, more Importantly, can be shown 
to be a transportation problem as was first pointed out 
by Bowman (Bowman 56). Specifically, let 

^1J“ production in period i for demand 
in period J, J - 1 
then P4 can be written as 



N N 

Pii' : min Z? = Z Z a^ -j 

4 1=1 j«i 
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3 • ^ 



f 







< + yi - 

+ yi - yi-i 



i ■ 2 



1-1 



» • • • » 






(Pil.2* ) 




[du + vj, 



J - 1 
J - 2 
J - N 



N-1 



X - 0 and INTEGER 



where 



Hj - Pi + J - 1+1 



N 



ail • Pi 



and X Is Integer by the properties of the transporta 
tlon problem. 2 



branch and bound algorithm emmeratlon would have to 
be done on only a subset of the problem variables, l.e. 
y. For each feasible y generated, corresponding to a 
node In the enumeration, a transportation problem can 
be solved. This, In effect, implicitly fathoms each time 
what would have been a significant portion of the 

^The lower summation range in (P4.1’) could be ex- 
tended to 1 and the upper range In (P4.2‘) to N If 
a^^j =00,VJ < 1, 1 = 2,...,N. 



The impact of this development is that. In a 



i 
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enumeration tree had branch and bound been tried on 
the entire problem. Furthermore, the number of trans- 
portation problems solved can be cut to a small fraction 
of the number of nodes actually visited by using suit- 
able tests and bounds. In effect, the branch and bound 
is a search for the minimum cost combination of a 
feasible y and its corresponding optimal production- 
inventory schedule. 

Lower bounds for y can be obtained directly 
from constraint set (P3«3) which cam be solved for y as 




(2.4) y2 - 2y^ + d^ - d2 - v^ 

- yi-2 ^i-i “ i = 3 ,...,n 

This set of lower bounds on y is integer since all 
terms on the right-hand side are integer. 

Upper bounds on the problem variables must 
exist since the model is bounded , The following lemma 
will give a set of global upper bounds on y. 

Lemma 2.1 Global upper bounds for y are 

(2.5) y, - J* 2 ci. - id. - (J-l)v_ J » 1,...,N 

J 1=1 ^ i«l ^ ° 




■m 



i 
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Proof: For J * 1 from (2.2) 



Yl - Wi - di + Vq + Wq 



and by inspection the right-hand side is maximized when 
takes on its maximum value. The maximum capacity in 
period 1 would be a level sufficient to produce all 
future demands in period 1. This value is 



and so 



or 



w 



1 



N 

I d 
i=l 



1 



w 



o 



max y, = ^ d. - d 

1 1=1 1 




N 

E di - 
i=l ^ 



^1 



In general let J = k, ke{2 ,3 , • • • >N } 

k 

yj^ = kW]^ + (k-l)w2 + ... + ^'^o 



Again the right-hand side is maximized when Wq^ is, and 
from the argument above 



< N k 

yj^ - k* E d4 - Id. 

^ i“l ^ i»l i 



(k-l)Vo 



so the result is true in general. 



QED 
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As will be discussed later In this chapter, 
upper bounds are not really crucial to the algorithm 
since more powerful termination rules will be employed. 
A special case which has an optimal solution with 
yj^ “ 0 will also be Investigated subsequently. 

11,3 Properties of P3 

Since a branch and bound algorithm Is to be 
developed It Is Important to have efficient bounds and 
termination rules In order to perform as little enumer- 
ation as possible. The following lemmas will establish 
useful properties of some variables of P3. 

Lenina 2.2 If all y^^ are fixed for 1 = 1 , . . . ,N-1 and y^^ Is 
Incremented by n to 

ne{l,2,3,. . 

then the new value of Is 

'll ■ "n ^ " 

Proof: From (2.3) 

" ^%-l ‘‘‘ ^N-2 ^ ~ ^-1 

= (y^+n) - 2yj^_^ + y ^_2 + ‘ Vl 

" QED 

For Lemmas 2 . 3 and 2 . ^4 , Corol lary 2.2.2 and Theorems 
2.1 and 2.2, when given a set of y^^, 1 = l,...,j, it Is 



^ 0 . 



assumed that the remaining are at their lower bounds for 
k = 

Lemma 2,3 If for given y^, 1 » the lower 

bounds of yj+j^ are positive, k = then * 0, 

k ■ 1, . , . ,N~J . 

Proof: For any ke{ 1,2, . . . ,N-J } from (2.4) and the 

hypothesis 



yj+k " 2yj+k_i - yj+k-2 «ij+k-l ' <^J+k ^ ° 



Further, if 
on the left 



yj+k lower bound 

The expression for Wj^j^ 



, equality holds 
from (2.3) is 



’'j+k “ yj+k - ^yj+k-1 + yj+k-2 + ^J+k-1 - dj+k 



but the right-hand side has Just been shown to be zero. 
Thus 

“j+k = ® 

QED 

Lemma 2 . 4 If for given y^, 1 * 1,...,J, the lower 
bounds of yj+k positive, k » 1,...,N-J, and y’ = 
yj + n, ne{ 1 ,2 ,3 , . . . } then 

w' * w. + n 

J J 



! 

I 



Is 






and for k * the lower bound on y' 

j+k 



yj+k + (k+l)n 



and 



w 



J+k 



0 



Proof: Prom (2.3) 



' yj - 2yj-i ^ yj-2 * - "“j-i 

SO 

Wj - (yj+n) - 2yj_;^ + yj_2 + - dj.^ - Wj + n 

which proves the first assertion. The second and third 
are proved by Induction on k. For k = 1, from (2.4) 



'j+i ■ ^ 



2n 



and from (2,3) and Lemma 2,3 



"j+i ° ■ 2(yj*n) + yj-1 + <ij+i - dj 



= = 0 



Assume the relationship holds for a general ke{2,3>...» 
N-J-1). Then for k+1 from (2,4) 



i 

I 
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'’j+k+1 ■ - [yj+k-l'^'^'’^ + dj+k - dj^.k+1 

“ yj + k+1 ■*■ (k+2)n 



The final assertion is then true by Lemma 2.3. 

QED 

Remark: From the relationships (2.3) and (2.^) it can 

be seen that a change in leaves and y^, 

..•»yj_l unaffected. 

Now let 

C(yj) “ capacity cost as a function of yj 
Pl(yj) “ production-inventory cost as a 
function of yj 

T(yj) * total cost as a function of yj 
so 

T(yj) = C(yj) + Pl(yj) 



Then the preceeding lemmas can be used to establish 
the following cost properties for P3: 

Theorem 2. 1 If for given y^, i = 1,...,J, the lower 
bounds for yj + j^ are positive k = 1,...,N-J then C(yj+n) 
is linear and non-decreasing in n for ne{l,2,3. . . }. 

Proof: From Lemma 2.H for y* = y.+n, ne{l,2,3...) 

J J 



w' “ w. + n 

J J 






k = 1,2,. . . ,N-J 



0 



Thus, from PI 

N 

C(y +1) - C(y.) “ C(y +2) - C(y. + 1) =... = (c.+ I m. ) - 0 
J J J J J 1=J 

QED 

Before stating the next theorem, a well 

known result for the linear program associated with 

P4' must be mentioned. Since the linear program 

satisfies the sufficient conditions for unimodularity, 

then all extreme points of the convex hull are integer 

(Hu 69 ). This allows use of the properties of the 

general linear program in the following proof. 

Theorem 2.2 If ne{l,2,3,...} and for any Je{l,2,...,N} 

if y' = y. + n and for each value of n the y. are set 

at their lower bound, 1 = J+1,...,N, then Pl(y') is 

J 

convex in n and the y^ , i = J+1,...,N. 

Proof: Assume the following linear programming 

problem 



(pi) : 

s . t. 



rain c*^x 



Ax = b 



X 



> 



0 



and the two parameterizations of b 



I 



1 

1 



( 
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(p2) 


Ax - b + (l-A)e 


(p3) 


Ax = b + e 



where 0 - A - l,e Is any vector which maintains feasi- 
bility in (p 3) and the optimal solution of (pi) is 
and of (p3) is where 

0 *^X 2 ^ = f(b) 

0 *^X 2 * f(b+e) 

are the optimal objective function values as functions 
of the right-hand sides. Now let 

(2.6) x^ = + (1 -A)x2 0 - a - 1 

and since the problems are linear it must be that 

(2.7) 0 *^X 2 “ 0*^2^ + (1- A) 0 *^X 2 

Multiplying the constraints of (pi) and (p3) by the 
appropriate constants gives 




(1-A)Ax2 ■ 



(1-A)(b+?) 
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and since X is a scalar, this can be rewritten as 



Adding these two constraint sets gives 

A[XXj^ + (l-X)^^] = b + (l-A)e 
which, from (2.6), is Just 

Ax^ - b + (l-X)e 

which proves that the convex combination of the op- 
timal solutions to (pi) and (p3) is feasible to the 
convex combination of right-hand sides. Now if x^ is 
the optimal solution to (p2) 

(2.8) c^x^ - c*^x^ 

since it is a minimization problem. 

From (2.7) and (2.8) 



AXx^ « Xb 



0 - X - 1 



A(1-X)x 2 ' (1-X)(b+S) 




Xc'^Q^ + (l-X)c'^X2 



0 - X - 1 
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producing the general convexity result and, by the 
remark preceding this theorem, the desired result for 



PI(yj+n) . 



QED 



The following corollaries follow trivially. 

Corollary 2.2.1 PI(yj^+n) is convex for ne{l,2,3, . . . } . 
Corollary 2.2.2 If for given y^^ , i * 1,...,J, the 
lower bounds of are positive, k ^ 1,...,N-J then 

PiCy.+n) is convex in n for ne{l,2,3,. . . 

Proof: From Lemma 2.4 the capacities in periods' 

J+1 through N are equal to the capacity in period J. 
The result follows from this and Theorem 2.2. QED 
The major result can now be given as follows: 

Theorem 2.3 If for given y^, 1 = 1,...,J, the lower 
bounds for yj+j^ are positive, k = 1,...,N-J, then 
T(y +n) is convex In n for ne{l,2,3, . . . } . 

Proof: The proof follows from Theorem 2.1 and 

Corollary 2.2.2 and the fact that the sum of convex 
functions is convex. qED 

Corollary 2.3.1 TCy^^+n) is convex for ne{ 1,2 ,3 , . . . } . 

II. 4 Optimality Tests and Bounds for the Enumeration 



The results of section II. 3 Ccin now be used 



to develop bounds and optimality tests for the 



^7 

algorithm. The logical first step is to solve the 
linear programming problem for PI and use 
Optimality Test If the optimal linear programming 

solution is all Integer, then this is an optimal solu- 
tion to P3. 

Otherwise let 

Z* - optimal objective function value of the 
LP 

linear program associated with the integer 
programming problem 

Z*p = optimal Integer program objective function 
value 

It is a well known fact that 



Z* ^ 
LP 




and this is used as 

Bound 1; A lower bound on total cost is Z* . 

LP 

The algorithm will use Bound 1 in the follow 
Ing optimality test. 

Optimality Test 2: If Zjp is the current integer pro- 

gramming solution and 

then Zjp Is an optimal solution to the Integer 
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program, 3 

The next bound uses the minimum production- 

inventory cost, PI*, to obtain an upper bound on 

capacity expansion cost. PI* can be obtained by 

solving P4' with the values of the right-hand sides 

of constraint set (P4.1’) equal to infinity. This 

corresponds to putting an infinite supply at the 

sources of the transportation problem. 

Lemma 2.8 If Z1 is a current optimal solution in the 
3 

enumeration then 

a = Z» - PI* 

3 

is an upper bound on capacity expansion cost. 

Proof: Assume there is a capacity expansion cost 

01^ such that 

> cx 

and with a corresponding such that 

3 

Z^ < Z' 

3 3 



But this implies that 



3j”a^ indicates the next integer greater than a. 



I 



PI 






al < 



Z' - 

3 



a = PI« 



^ 9 . 



which is a contradiction of the construction of PI*. 

QED 

Using the above gives: 

Bound 2^: If Z^ is a current optimal solution of the 

integer program, then an upper bound on capacity ex- 
pansion cost a is 



a - Z' - PI* = a 
3 



Lemma 2.8 can be useful in another way. If 
the current upper bound on capacity cost is a then 
certainly an upper bound on period 1 capacity is^ 



< 

”l * 



a 



N 

Ci+ Z m^ 
J- i*i 1 



since any larger value of Wj^ would violate a, A 
similar bound for W2 is 



wp - 



- / N 
a- (ci + Z m. )w, 
-L i=i 1 1 

N 

Cc2+^£2^i^ 



4 




indicates the greatest Integer smaller than a. 
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In general the following bound can be used. 
Bound 2* the current a, upper bounds on capacity 
expansion In each period are 




/V 



N 



ct + Z m^ 

1=1 i 




J “ 2, 



,N 



The following bounds use the results of the 
previous sections. 

Bound For any feasible values of y^^, 1 » 

lower bounds on the values of y^^, 1 » J+1,...,N can be 

obtained from (2.4). 

Bound 5: Global upper bounds on y^, 1 « 1,...,N can 

be obtained from (2.5). 

Optimality Test 3 : If enumeration has been completed 

with y^^ at Its global upper bound or the upper bound 
Inferred from Bound 3» whichever Is lower, then the 
current optimal solution Is optimal for the Integer 
program. 

The following bound on yj^ uses the results 



of Corollary 2.3.1. 
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Bound If for a given set of 1 = 

T(y„) - T(y^+l) 

then the enumeration on for this set of yj^ can be 
fathomed at (y^^+i;. 

The following bound uses the result of 
Theorem J.3. 

Bound 7 : If for a given y^, lower 

bounds on y^ are positive for 1 = J+1,...,N and If 

T(yj) - T(yj+1) 

then the enumeration on yj can be fathomed for this set 
of y^^, 1 = l,...,j-l at (yj + lj. 

Bounds 8 and 9 make use of Theorem 2.1 and 
Corollary 2.1.1. 

Bound 8: For given y^^, 1 = l,...,j. If the lower bounds 
on are positive for 1 = j+l,.,.,N and If C(yj) - a, 
then yj^ Is fathomed for the preceding set of y^. 

Bound 9.: If for given yj^, 1 = 1,...,N-1, C(yj^) - a, 

then no further enumeration on yj^ Is necessary. 

Since the algorithm conducts the search by 
enumeration on yj^ most often. Bounds 6 and 9 become very 
powerful tools. Bounds 7 and b are also extremely 
powerful because their conditions will be satisfied 
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often during enumeration on the later period variables 
due to relationships (2.4). The usefulness of these 
bounds diminishes the need for the global upper bounds 
(2.5) since such bounds must include all cases where 
Bounds 6 through 9 would be useful. 

II. 5 Algorithm for P3 

Step 0^ Initialization Set the value of current optimal 



cost, Z^, and the current bound on capacity cost, a, to 
infinity . 

Step ^ Solve the associated linear program. If the 
solution is integer, stop. (Optimality Test 1) Other- 
wise set the linear programming bound to (Bound 1) 



and go to Step 2a. 

Step 2a Solve an uncapacitated transportation problem 
to get the minimum production-inventory cost PI* and go 
to 3tep 2b. 

Step 2b Find global upper bounds for y^, i = l,...,N 
from (2.5). Oo to Step 3. 

Step ^ Start the enumeration by setting 
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= max [ ( w^+v^-d^ ) , 0] 

and finding lower bounds for the remaining values of 
i » 2,...,N, from (2.4), using this value of y^^ . The 
pointer (PNT) is N. Go to Step 5. 

Step 4a PNT = PNT - 1. Go to Step 4b. 

Step 4b With PNT = J, let yj = and, if J < N, find 

new lower bounds for succeeding y^^ from (2.4). If J 1 
and any lower bound exceeds an upper bound (from Bounds 
3 and 5), yj is fathomed so go to Step 4a with PNT = J. 
If 1 and some lower bound exceeds an upper bound, go 
to Step 10 (Optimality Test 3). Otherwise go to Step 5 
with PNT = J. 

Step 5 With PNT = J determine capacity cost a = C(yj). 
If a > 3 and either J * N, or J < N and yj^ > 0, 1 » J+1, 
...,N go to Step 4a v/ith PNT = J (Bounds 8 and 9). 
Otherwise if a > 3 go to Step 4b with PNT = N. If 
0-3 and if J < N with > 0, i = J+1,...,N go to 
Step 6 with PNT * J (This sets up convexity test for 
Bound 7). Otherwise go to Step 6 with PNT = N (This 
sets up convexity test for Bound 6). 

Step 6 With PNT = J solve a transportation problem 
to obtain the related optimal production-inventory cost 
Pl(yj). Then total cost is T(yj) = PI(Yj). Go 

to Step 7 with PNT = J . 



5 ^. 



Step 7. With PNT *J, IfJ = NorJ <N and > 0, 

1 = and T(yj) > T(yj-l) then yj Is fathomed 

so go to Step 4a with PNT = J (Bounds 6 and 7). If 

T(yj) = T(yj-l) » In any case go to Step 9. For all 

other cases If T(y^) < Z* go to Step 8, otherwise go to 

3 



Step 4b with PNT = N. 
Step 8 If T(yj) * |~Z 



Z» 

LP 



go to Step 10 (Optimality 



Test 2). Otherwise go to Step 9* 

Step 9 Either Include T(yj) = Z^ on the list of 
current optimal solutions or replace the current list 
with Z^ = T(yj), whichever Is appropriate. Update In- 
termediate bounds from Bound 3 and go to Step 4b with 



PNT = N. 



Step 10 Decode the current optimal solutions and STOP. 

The algorithm, as written, will also find 
all alternative optimal solutions If termination Is not 
due to Step 8. The flexibility of choice among solu- 
tions would normally be desirable In a management 
situation. If the alternate optima and longer running 
times are not desired then the following changes to the 
algorithm are needed: 

Step 3 With PNT = J determine capacity cost a = C(yj), 
If a - S and either J ** N or J < N and y^^ > 0, 1 = J+1, 
...,N, go to Step 4a with PNT = J. Otherwise If a - 3 
go to Step 4b with PNT = N. If a < a and If J < N with 
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> 0, 1 = go to Step 6 with PNT = J, 

Otherwise go to Step 6 with PNT * N. 

Step 7 With PNT » J, If J = N or J < N and > 0, 
i = J+1,...,N and T(yj) - T(yj-l) then yj Is fathomed 
so go to Step ila with PNT *= J . If T(y^) < Z’ go to 
Step 8, otherwise go to Step 4b with PNT = N. 

Step 9 Replace the current optimal solution with 
= T(yj). Go to Step 4b with PNT » N. 

II. 6 The CPI Model Under Special Demand and Cost 
Assumptions 

In this section two special cases of the CPI 
model will be discussed, one Involving monotone non- 
decreasing demands and the other monotone non-lncreas- 
Ing demands. The former leads to significant simpli- 
fications In the CPI algorithm, while the latter Is a 
trivial case of the CPI model. 

The following two cost assumptions hold In 
this section: 

(a) The cost assumptions of PI hold. 

(b) The cost coefficients of PI (Chapter I) are mono- 
tone non-lncreaslng with time. 

Assumption (b) allows the possibility of discounting. 
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Monotone non-decreasing demand 



The case of demand which is monotone non- 
decreasing is likely to occur often in the context of 
firms Interested in solving capacity expansion problems. 
The main result for this case is contained in the 
following theorem. 

Theorem 2 . 4 If assumptions (a) and (b) hold, demand 
is monotone non-decreasing, and 



Nw + 
o o 



< 



N 

I d 
i=l 



1 



then there is an optimal solution with Yfj “ 0. 
Proof: Assume in contradiction that 



= a > 0 

It can also be assumed, as discussed in Chapter I, that 
Vj^ * 0. But this Just implies that there are a > 0 
idle units of capacity in the model. Consider the first 
unit of idle capacity. One of two cases can be assumed 
to exist. 

Case I- The idle unit was added in period J instead of 
k where it was first used, where J < k ^ N. 

Case II- The capacity was added in some period previous 
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to .5, was used up to and including J-1, and became 
idle in J . 

For Case I to happen either the cost of adding capacity 
in J must be cheaper than k, i.e. 



N 



N 



c^ + E m^ < Cw + Em. 

= . i k i 



•J iSj 



which, since (assumption (b)) implies that 

k-1 ^ 

^E m^^ < 0 contradicting assumption (a); or the cost is 
equal to that in k in which case the idle unit could be 
removed in J and added in k with no change in cost. 

For Case II it must be that Wj * 0 since if not, then 
Case I applies. So this unit of capacity must have 
been built in some previous period q and used to pro- 
duce from q through J-1. Now at least one of these 
(J-q) units of production must be inventoried into 
period J since if Vj_]^ = 0 then since > Xj this 

implies that ^ dj contradicting the demand assump- 

tion. Furthermore, by the same argument, v^ - 1, i = 
q,q+l, . . . , J-1 so the idle unit in J was used to produce 
for Inventory in period q and at least that unit of in- 
ventory was carried into period J . The production 
schedule of periods q+1 to J-1 will not change if the 
unit of capacity is added in q+1 instead of q, however 
the unit which was Inventoried from q to J could be 
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produced in J instead, eliminating the idle unit of 
capacity. With idle capacity the Incremental cost con- 
sidering .lust the one unit of production in question 
was 



N J-1 

(2.9) c + E mi + p + Eh. 

q i*=q 1 q i=q i 

while adding the capacity in q+1 and producing in J 
gives 

(2.10) Cq+j^ + + Pj 

Subtracting (2.10) from (2.9) gives 



(Cq-Cq+i) 



J=l. 



+ "»q + ^Pq-Pj ^ - 0 



by assumptions (a) and (b). Thus, in both cases, idle 
capacity can be reduced with, at worst, no reduction in 
cost . 



If idle capacity still exists, the argument 
can be repeated. 



QED 

Corollary 2.4.1 ; If the assumptions of Theorem 2.4 
hold then all capacity is used in each period. 

The use of Corollary 2.4.1 leads to an Imme 
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dlate simplification of the production-inventory prob- 
lem . Since the right-hand sides of constraint set 
(P4.1) are Just capacities, equality holds in all N 
constraints and the are determined trivially for a 
given y. Constraint set (P^l. ki) becomes N equations in 
N unknowns (since v^ is explicit and Vj^ is either zero 
implicitly or set explicitly) and is thus deterministic. 
So the subproblem in the enumeration can be solved by 
inspection each time. 

Another observation follows from Theorem 2.4 
and Lemma 2.4.1. If all capacity is always used it can 
certainly never be feasible to build capacity in any 
period sufficient to produce all remaining demand. In 
fact, the upper bounds on capacity Increase are much 
smaller than those given by (2.5). Instead of stating 
them explicitly, they can be handled implicitly in the 
algorithm by the following: 

Bound 10 ; For given y^^, i = 1,...,J, if the increment 

yj = yj + 1 

causes the succeeding lower bounds to contain yj^ > 0, 
then yj is fathomed for the preceding set of yj^. 

Bounds 6 through 9 are now useless, since 
their conditions are never met in this case, so 



I 

I 
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Bound 10 takes over as a main fathoming rule In the 
algorithm. Furthermore, the global upper bounds from 
(2.5) are Irrelevant to this special case and can be 
Ignored . 

The basic algorithm requires the following 
modifications; 

1) Drop Step 2b. 

2) Step 4b* With PNT = J let yj = yj+1 and find new 
lower bounds for the succeeding elements from (2.4). 

If yj^ > 0 then yj Is fathomed so go to Step 4a with 

PNT = J (Bound 10) . If J =1 and y^ exceeds Its upper 

bound from Bound 3 or y^j > 0, go to Step 10. Otherwise 

go to Step 5* with PNT = J. 

3) Step 5 * With PNT = J determine capacity cost a = 
C(yj). If o > go to Step 4b* with PNT = N-1. Other- 
wise go to Step 6* with PNT ® J. 

4) Step 6 * With PNT = J, determine the related optimal 
production-inventory schedule by Inspection and the re- 
lated cost Pl(yj). Total cost Is T(yj) = C(yj) + Pl(yj). 
Go to Step 7* with PNT = J. 

5) Step 7* With PNT = J If T(y. ) > Z* go to Step 4b* 

J 3 

with PNT = N-1. If T(yj ) « go to Step 9*. Other- 
wise go to Step 8. 

6) Step 9* Either Include T(yj) = on the list of 
current optimal solutions or replace the current list 
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with = T(yj), whichever Is appropriate. Go to Step 
4b’ with PNT = N-1. 

Monotone non - Increasing demand 

Assume that demand Is monotone non-lncreaslng 
and that monotonlclty Is not destroyed by Initial In- 
ventory, l.e. 



d 



1 



- V, 



> 



d 



2 



By assumption (b) production should always take place 
as late as possible. Then for the CPI model the only 
optimal capacity expansion policy can be to expand 
capacity In the first period to satisfy the demand In 
that period, l.e. 



«1 = <*1 - Vo - 

This will then be sufficient to produce the remaining 
demands, so no capacity expansion problem need be 
solved. If the cost assumptions do not hold however, 
then this case Is Just handled with the general CPI 



model 
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II. 7 Transformat Ion for P2 



A similar procedure to that used In trans- 
forming the CPI model PI will be used to transform the 
CIDPI model P2 . First a lemma will be proved which 
will be useful In Interpreting the transformed model, 
and In the algorithm to follow. 

Lemma 2 . 9 Por any period 1 either 

(1) w^ > 0 and r^^ = 0 
or 

(2) Wj^ = 0 and r^^ = 0 
or 

(3) w^ = 0 and r^ > 0 

Proof: The proof Is by contradiction and con- 

struction. Assume that w^ > 0 and r^^ > 0 for any peri- 
od 1 gives a minimum cost. There are three cases: 

Case 1 w^ - r^ > 0 

Let w' = w^ - r. and r' = 0 

1 ^ i 1 

Then 



c j^w^-f j^r |+( w^-r^) 



N 

Em. 

J = 1 J 



N 

Ci(Wi-ri)-fi* 0+(wi-ri) -^E^mj 



N 



:iWi- Cj^ri+(wi-ri)-^E^mj 



N 



= 1^1- + 



< 
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since c^> 

The assumption Is contradicted and (1) has 

been constructed. 

Case ^ - r^ = 0 

Let w’ =0 and r' = 0 
1 1 

Then certainly 
0 < 



since c^ > 

The assumption Is again contradicted and (2) 

has been constructed. 

Case ^ - rjL < 0 

Let w’ = 0 and r’ = r. - W4 
1 111 

Then 



N 



c ji^r ' + (wj^-r ^ nij = c j^(r j^-Wj^) + (wjL-r • 

1 "" 1 J * 



N 



flWi-firi+(wi-ri)‘ l^mj 



N 






since CjL > 

The contradiction Is complete and (3) has 
been constructed. 



QED 



2SW U 



I 
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Now using Standard Transformation 1 as in 
Section II. 2, constraint sets (P2.1) and (P2.2) can be 
combined to give the set of relations 



(w2-r2) = d2 - di + 72 - 2y^ + Vq 

( 2 . 11 ) 

(w3-r3) = d^ - d2 + 73 -272 + 7 i 
(wj^-r^j) » - dj^_i + 7 n - 27N-1 + yN-2 



Note here that, b 7 Lemma 2.9 and w, r - 0, If the right- 
hand side of a constraint in ( 2 . 11 ) Is positive, then 
> 0 and r^ * 0. If it is zero, then Wj^ = r^ = 0, 

If It is negative, then = 0 and r^ > 0. 

Solving (2.11) for r and substituting In (P2) 

gives 



N 

P5: min Z 5 = ^^^CPlXi + h^Vi + f j 7 i + (ci-fj^)wj^] + M 



s .t . 



xi ^ di+ 7i - Vq 

(P5.1) . 

Xi - di + 7i - 7i-l i = 2,...,N 



n 






I 



1 

I 

i 

i 

I 



I 

( 



i 
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(P 5 . 2 ) 



(P 5 . 3 ) 



(P 5 . 4 ) 



Xi - vi = di - Vq 

Xi + =d^ i = 2,...,N-1 

""N + '"n-I “ 

di + yi - Vq - 0 

di + Yi - Yi-i - 0 i = 2,...,N 

wi - di + Yi - (Wq+v^) 

W2 - d2 - d^ + Y2 - 2 y3^ + 

- di - dj^_3^ + Yi - 2yi_i + yi_2 1 = 3,..,N 



X, 7 , W, y ^ 0 and INTEGER 



where 



- 2 fj ^^2 ■*■ ^ 1+2 " ^i +1 ^ “ l,...,N -2 



^N-1 “ 






M 



= V(f^-fi^^)d^ + fj^d^ + - f^Wo + (f2-fi-mi)v^ 



Constraint sets (P 5 . 1 ) and (P 5 . 2 ) and the 
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first two sets of terms In the objective function are 
Just problem P4 and thus are equivalent to a transporta 
tion problem. Writing constraint set (P5.^) separately 
with the last set of terms of the objective function 
gives 

N 

P6: min Zg = ^^^(Ci-fi) Wj^ 



s . t . 

wi - - Wq 

Wn - - d^.i + Yn - 2yj^_i + y^ -2 

w ^ 0 and INTEGER 

which, given y, is solvable for w by inspection since 
the coefficients in the objective function are positive 
and the right-hand sides of the constraints are integer. 
Furthermore, it is possible that the right-hand side of 
some constraint i is negative or zero. Recalling Lemma 
2.9, it follows that in this case “ 0 and r^^ equals 
the absolute value of the right-hand side. So, for 
each possible y in the CIDPI model, it is necessary to 
consider a transportation problem and a trivial linear 



program 
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Essentially the same branch and bound enumer- 
ation scheme can be followed here as was used In the 
capacity Increase case. Lemma 2.1 still applies as a 
finite bound on the enumeration and (2.5) Is still a 
valid set of global upper bounds. 

The lower bounds on y are of a slightly dif- 
ferent form than (2.4) and can be gotten from (P5.3) and 
7 - 0 as 

yi - Vo - di 

(2.12) - yi_i - dj^ 1 = 2,...,N 

7 > 0 



II. 8 Pr'opertles and Algorithm for P5 

The properties of P5 are similar to those of P3. 
Because of the new lower bounds (2.12) however, most of the re- 
sults must be restated and reproved. 

The statement and proof of Lemma 2.2 hold equiva- 
lently for P5 if the change of variables 

= (w^ - r^) 1=1,... ,N 

Is made. To arrive at sufficient conditions for convexity to 
hold In the optimal values of Z^, the remaining results must be 
reiterated. For Lemmas 2.10, 2.11 and 2.12 and Theorems 2.5 and 



I 
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and 2.6, when given a set of y^, 1 = It is assumed 

that the remaining are at their lower bounds for k»l,...,N-J. 
Lemma 2.3 is restated as 

Lemma 2 . 10 If for given y^^, i = 1,...,J the lower 
bounds of yj+k sire positive, k * 1,..., N-J then 
* 0 f'or ^ ® 2,..., N-J. 

Proof: For any ke{ 2, 3, . . . ,N-J} from (2.12) 

yj+k " yj+k-1 - «ij+k > 0 

and 

yj + k-1 “ yj+k-2 - ^J+k-1 > 0 

where equality holds on the left. Then from (2.11) 

('^J+k-^J+k) = (yj+k-yj+k+l+dj+k)-(yj+k-i"yj+k-2'*‘^j+k-l) 



Next, Lemma 2.4 becomes 

Lemma 2.11 If the conditions of Lemma 2.10 hold and 
y’ ■ yi+>^ net 1, 2, 3, • • • ) then 

(Wj-rj)' = (wj-rj) + n 

and for k = 1,2,..., N-J the lower bound on y' 

J+k 



is 
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yj+k + ^ 

and 

and 

('^j+k-rj+k)’ = 0 k = 2,...,N-J 

Proof: Prom (2.11) 



(Mj-rj)’ = dj - dj.i + (yj+n) - 2yj_i + yj.j 
= (wj-rj) + n 

For k = 1 



^j + 1 “ ■* " (yj+n) - dj+i = yj+i + n 

so from (2,11) 

^'^j+l~^J+l) ' = ^J+1 ” (yj+l+n) - 2(yj+n) + yj_^ 

= (wj+i-rj+i) - n 



Doing induction on k, for k = 2 

y * ~ y ' “ 1 4-0 

J+2 J+1 J+2 



= (yj+i+n) - d 



J+2 



= yj +2 + 
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Assume the relationship holds for a general k. 

Then for k+1 



^J+k+1 “ ^J+k ” ^J+k+1 

= (yj+k+n) . 

= ^J+k+1 ^ ^ 

Finally by Lemma 2.10 

^'^J+k”^J+k^ ' ^ ° k = 2,...,N 

QED 



The following Is an additional property which will be 
needed. 

Lemma 2.12 If the conditions of Lemma 2.10 hold, then 
Proof: Prom (2.11) 






since y' Is at Its lower bound so the first term Is 
zero, and the second term must be non-negative by (2.12) 

QED 



Now Theorem 2.1 can be restated as 
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Theorem 2,5 If for given 1 = the lower 

bounds on are positive for k = 1,...,N-J, then 

C(y.+n) Is strictly convex In n for ne{l,2,3> • • • ) • 

J 

Proof: For yj = ^>3, . . . >from Lemmas 

2.11 and 2.12 



(wj-rj ) ’ « ” 

^"^J+l-^J+1^’ ' (wj+l-^j+l) - n ^ 0 

+k“^j+lc^ * * 0 k * 2,...,N 

The capacity change variables In period j could start 
from rj > 0 and go to Wj > 0 as n Increases. Assume 
there Is an n^^ such that r^ > 0 for n < n^^ and Wj > 0 
for n > n^^. Then the change In capacity cost for each 
Increment of n < n^^ Is 

Ai - fj + mj - fj+i 

while for n - n^^ the cost change Is 

. Cj + mj - > 0 

where A 2 ^ 0 follows from the cost assumption of P2 
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that 

""l Ji"*k + Jj^k VJ > 1. 1 - 1,...,N-1 

Further from P2 it was assumed that Cj > fj giving the 
convexity result 



Corollary 2.5.1 C(y^+n) is convex in n for ne{l,2, 

3, . . . } . 

Theorem 2.2 and Corollary 2.2.1 still hold 
since, as has been pointed out, the production-inven- 
tory section of the model is equivalent. The main 
result. Theorem 2.3* must be restated as 
Theorem 2 . 6 If the conditions of Theorem 2.5 hold, 
then T(yj+n) is convex in n for ne{l,2,3, . . . } . 

Proof: The proof follows from Theorem 2.5 and 

Corollary 2.2.2 and the fact that the sum of convex 
functions is convex. 

QED 

Corollary 2.3.1 is unchanged. The above properties of 
the P5 variables and convexity in the optimal values of 
will, with a few exceptions, permit the bounds and 
optimality tests of section II. 4 to be used. 
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Algorithm for the CIDPI model 

The bounds and optimality tests of section 
II. ^ hold for the CIDPI model with a few exceptions. 
Bound 3 Is no longer valid since the possibility of 
decreasing capacity allows a w^ to exceed these bounds 
and not violate a. For Bound h the lower bounds are 
now found from (2.12). Optimality Test 3 can no longer 
use Bound 3. Finally, since Theorem 2.5 only gives the 
weaker result of convexity for C(y), Bound 8 has to be 
rewritten to ensure that the test Is being conducted 
in the area where the slope of C(y) Is positive. 

Bound 8 * For given y^^, 1 = l,...,j If y^^ > 0,1 = J+1, 
C(yj) - C(yj+1), and C(yj+1) - a then yj is 
fathomed for the preceding set of yj^. 

The basic algorithm of section II. ^ requires 
some changes for the CIDPI model. The lower bounds In 
Steps 3 and 4b are now found from (2.12) and the Initial 
value for y^ In Step 3 Is 

y^ = max [(vQ-djL^* 

The first part of Step 5 must be rewritten to reflect 

Bound 8* as follows: With PNT = J determine capacity 

cost a = C(y.). If a > a and either J = N or J > N, 

J 




I 



I 



€ 
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> o for 1 “ with C(yj) - C(yj-l) go to 

Step 4a with PNT = J (Bounds 8’ and 9) The remain- 

der of Step 5 Is the same. Finally, the update of 
Bound 3 In Step 9 can be dropped. 

The remarks of section II. 5 pertaining to 
alternate optimal solutions are relevant here also. 
Finally, note that lower bounds (2.12) will serve to 
Increase the feasible space for y and thus longer 
average running times should result for the CIDPI model. 

II. 9 Cases of Monotone Demands with Special Cost 
Structures for the CIDPI Model 

As with the CPI model, the efficiency of the 
algorithm for the CIDPI model Is Increased greatly when 
monotonlclty exists In demand. The following hold for 
this section 

a) The cost coefficients of P2 (Chapter I) 
are monotone non- Increasing with time. 

b) The cost assumptions of P2 hold. 

Monotone non-decreasing demand case 

Lemma 2.13 If assumptions (a) and (b) of this section 
hold, demand Is monotone non-decreasing and 



I 




I 
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Nw + V - 
o o 



N 

I d 
1=1 



1 



then there is an optimal solution with yjj = 0. 

Proof: The proof follows trivially from the fact 

that these are the same conditions as for Theorem 2,^4. 
Now, In addition to the possibility of eliminating 
idle capacity by the arguments of Theorem 2.^, there is 
the added possibility of scrapping it and recovering a 
salvage value, so the results hold a fortiori, QED 

The conclusions of section II. 6 which follow 
from Theorem 2 A are valid here for the result above; 
i.e. the transportation subproblem becomes trivial and 
can be solved by inspection for a given y, and a much 
more efficient set of upper bounds can be included 
implicitly in the algorithm. That is. Bound 10 of sec- 
tion II. 6 can be used instead of Bounds 6, 7> 8, and 9. 
The algorithm changes of section II. 6 remain valid with 
the exception that the new lower bounds in Step ^4b’ are 
found from (2.12). 



Monotone non-increasing demand 

For this case it must be assumed that 




in order that the monotonicity be preserved. 



I 
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Lemma 2.1^ If demands are monotone non-lncreaslng, 
assumptions (a) and (b) hold and further 

k-1 

(c) £ "'l > J> J = 

^ J 

then = dj^ - Wq - and = 0, i = 2,...,N. 

Proof: The statement for results from the fact 

that no backordering Is allowed and by cost assumptions 
(a) and (b) and monotone demands. It can never pay to 
build extra capacity to produce early. Assumption (c) 
eliminates the possibility of scrapping capacity In any 
period to the extent that capacity must be rebuilt 
later, which, by monotonlclty of demand, produces the 
result that w^ = 0, 1 = 2,...,N. 

QED 

Lemma 2.15 If demands are monotone non-lncreaslng and 
f - 0 then y^^ = 0. 

Proof: The proof follows trivially from the fact 

that demand Is monotone non-lncreaslng so, with f - 0, 
any Idle capacity In a period can be salvaged for 
further revenues. 

QED 



As was pointed out In section II. 6, when 
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= 0, all capacity is used. The transportation sub- 
problem of P5 thus can be solved by Inspection for a 
given y, speeding up the CIDPI algorithm greatly. The 
remaining remarks on algorithm changes for this case 
are the same as for the preceding case of monotone non- 
decreasing demands. 



Chapter III 



In this chapter, a special case of the CIDPI 
model will be presented which leads to a formulation of 
the entire model as a network. Other extensions of the 
basic model are then presented to allow piecewlse- 
linear, convex production costs and discontinuities in 
capacity acquisition. 

III.l Network formulation of the CIDPI model 

The formulations of Chapter II required a 
problem decomposition approach to provide suitable al- 
gorithms which took advantage of the problem structure. 
There is a special case of the CIDPI model, however, 
where the formulation can be simplified to that of a 
network and solved with a transshipment algorithm. 

This occurs when the cost of capacity increase is equal 
to the salvage value in any period. 

Such a model might be useful when there is a 
large second-hand market in the units of capacity. A 
firm might then acquire or sell units with equal ease in 
any period at the going market price. Maintenance 
costs then reflect the investment needed to keep each 
unit at its rated capacity for possible future sale. 
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The derivation of this case comes from recon- 
sideration of the CPI model PI. If w - 0 Is dropped, 
allowing 15 to be unrestricted, the problem would still 
be bounded If capacity costs Cj^, 1 = 1,...,N were mono- 
tone non-lncreaslng, since from constraint (Pl.l) and 
X - 0 It must still be true that 

1 > 

jEqWj — 0 1*=1,,,,,N 

This implies that each w^ Is bounded below by the total 
capacity built up to period 1, l.e. 

> 1-1 

w^ - - L W4 1 = 1, . . . ,N 

J = 0 '' 

since If the costs are monotone, building and subsequent 
scrapping of capacity can never result In a profit. 

That Is, with c^j^ - 0 and monotone 1 = 1,...,N, the 
right-hand side will always be as small In absolute 
value as possible. In fact, as a limit, it could never 
pay to build more capacity In any period than what Is 
needed to produce all remaining demand In that period, 
since the excess capacity would be Idle and can never 
reduce the total cost. Thus, with the assumption, which 
will hold for the remainder of this section, that 
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a set of lower bounds on v7 can be written as 



> 



Wl - -Wo 



> ^ 

W2 - Vo - J^di 

w. i - E d. 






Now with unrestricted in PI, if Standard 
Transformation 1 Is performed, then P3 no longer con- 
tains constraint set (P3.3). The Integrality of y 
still holds, however. To transform the remainder of P3 
Into a network formulation, slack variables s - 0 must 
be added to constraint set (P3.1). Then the sum of set 
(P3.1) Is subtracted from the sum of set (P3.2) to give 
a redundant constraint. Finally, multiplying the con- 
straints of set (P3.2) by -1, the transformed problem 
can be written as 



P6; 



min Zg 



N 



3. t. 

(P6.1) - yj^ + S3^ - dj^ - v^ 

Xj^ + ” ^1 ^ ®1 " ^1 ^ ~ 2 , . . . ,N 
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-X, + Vt = -d, + V 

11 1 o 

(P6.2) -x^ - Vi_i + 1 = 2,...,N-1 

-^n - ^N-1 = -^N 

(P6.3) - J^si = 0 

y, s - 0 and INTEGER 

where the constraint (P6.3) is redundant, the and K 
are as given in P3» and V|^ has been incorporated ex- 
plicitly in dfj. 

The coefficient matrix associated with P 6 
satisfies the sufficiency conditions for total uni- 
modularity as described by Hu (Hu 69 ), This guarantees 
that, since the right-hand sides of the constraints are 
integer, the linear program associated with P6 without 
the redundant constraint (P 6 . 3 ) will terminate integer. 
P6 posses more useful properties however, since the 
problem as written can be represented by the network 
displayed in Figure III-l. Thus the problem can be 
solved using a transshipment algorithm to obtain the 
minimum cost Plow. 

The only evident problems with the network 
are the cycles involving y and s which, if negative in 
cost, cause unboundedness. However, since flow Sj^ is 
costless, the monotonicity assumption on the cj^ is 
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^The figures on the source and sink arcs represent 
the arc capacities. For the remaining arcs the figures 
represent the flow variables associated with the arcs. 



V 
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sufficient to prevent this as seen from the following 
lenuna . 

Lemma 3 « 1 If capacity costs c^ are monotone non-in- 
creasing for 1 » 1,.,.,N, then 

N > 

lij'i - 0 J - 

Proof; It is obvious that 

c^ + mj^ - 0 

For J = N-1 

^CN-l"2cN+mj^_2~"^^ (cjj+m]^) = (Cj^_^-Cjj) + - 0 

and by Induction for j = k, 1 - k < N-1 

l=k^l * ^^k~^k+l^ + mj^. - 0 

QED 

It Is obvious from Figure III-l that the 
only cycles are of the form 

yj\j* ^N— 1* •••* ~ lj«««>N 

which have non-negative costs by Lemma 3.1. This con- 
firms the fact that, since P6 Is obtained by a linear 
transformation on the bounded problem PI, it too must 
be bounded. Emphasis must also be given to the fact 
that the capacity change costs c|, incorporated In P6, 
are completely non-stationary and thus admit to dls- 
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counting or any other aberrant structure. 

Networks and the general CPI and CIDPI models 

The existence of unimodularity in the special 
case of P6 permitted use of a network algorithm to ob- 
tain an optimal Integer solution. The general models, 
PI and P2, could also be formulated as networks but 
equivalence of the networks to the original models re- 
quires sets of constraints external to the networks. 
Moreover, these constraints violate the unimodularity 
conditions, as shown by the following examples, and 
thus exclude the use of any existing network algorithm 
to gain an integer solution* These examples are both 
linear programs having non-integer solutions with non- 
integer objective function values (eliminating alterna- 
tive optima). 

Example 1 : Example for the CPI model 

For N»6, c » (20,20,19,19,18,18), p » (5, 5, 4, 4, 3, 
3), h » (1,1, 1,1, 1,1), IR » (3, 3, 3, 3, 2, 2), d - (3,7,9 , 
13,18,21), Wq - 1, Vo = 1 

LP optimum 

Xi <= 2, X2 = 7, x^ * X]| = X5 = xg = 15.25 
v^ = 6.25, v^ = 8 , 5 , v^ = 

'^2 5, W3 = 8.25 



5.75 
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Objective Function = 735.25 

Example 2^: Example for the CIDPI model 

The problem Is as in Example 1 with the exception 
of f » (8, 8, 7, 7, 6, 6), d = (3,5.8,111,9,7) 

LP optimum 

Xj^ = X2 = * XI] » X5 = 7.6, X5 » 7.0 

=■ 5.6, V2 " 8.2, = 7.8, VI] = 1.4, v^ = 0 

= 6.6 r 5 = .6 

Objective Function = 457.2 

III. 2 Block Capacity 

In this section, CPI and CIDPI models will be 
extended to handle the case of block capacity. This 
terminology embraces any situation where capacity can 
only be acquired In units which are each capable of 
producing large numbers of units of demand in each 
period, as opposed to the Implicit assumption In effect 
so far, that each unit of capacity produced a single 
demand unit per period. Block capacity, then, corres- 
ponds to the economists' problem of "lumpiness" In 
acquisition which typically causes serious difficulty 
In a marginal analysis approach. 

In the remainder of this section It Is assumed 
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that In any period 

1 Capacity Unit » k Demand Units 
where k is a positive integer and is stationary, and 
that initial capacity, w^, is expressed in units of 
capacity. ^ 

CPL-Block capacity model 

Working first on the CPI model, in terms of 
block capacity, (2.1) can be rewritten as 



Following the same development as in Chapter II, and 
using Standard Transformation 1 to solve for w gives 



kw 




k(Nwj^+. . .+Wj^) " - dj^ + 



+ dj^ - (Nkw^+v^) 





N ” ^-1^ 



^Depending on the source of the demand, management 



has the option of choosing the most appropriate unit for 
demand, e.g. box, pallet-load, car-load, etc. The 



87 . 



and the expressions for can now be substituted back 
into PI to get the transformed problem. In this case, 
substitution in the ^5 - 0 constraints gives somewhat 
expanded results. The constraint set corresponding to 
(P3.3) now looks like 

- CYi + d^ - - 0 and INTEGER 

Ic 

« 

I ■ ^^N-i ■*■ % - ‘^N-l^ ’ ° INTEGER 

Taking the first constraint and ignoring the integer 
requirement gives 

yi - “"o * ''o - ‘^1 
while the Integer requirement gives 

i [y^ t - (kw^tv^)] . 

where n^ is a non-negative integer. But this is the 
same as 

y-=n.k-dT+kw +v 
11 1 o o 



value of k can then be set appropriately. 



I 

i 






A similar development for period J gives 
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- 2yj-l - yj-2 * dj.i - dj 

and 

- "'J-2 * ‘^J-1 - 

for n^ - 0 and integer. So, in general, substitution 
in w - 0 gives two sets of constraints, the first sim- 
ilar to set (P3.3) 



(3.1) 



y - kw + 
1 o 



V 

o 



d 



1 



i 2y^ * <^ 1 - <^2 ‘''o 
^1 - 2yj._i - yi.2 + di.i 



y 



> 



0 



1 » 3. 



,N 



and the second, a set of integrality constraints on 



y, = n,k - d- + kw + v 
11 1 o o 



( 3 . 2 ) 22 1 



n k + 2y + d, - d^ - v 



* ^^N-1 - ^N-2 ^N-1 



n^ i 0 and INTEGER 

1 » 1,...,N 



The effect of k, then is to superimpose a lattice of 



<< I 



interval k on the space of y. Up until now, k » 1 and 
30 the enumeration on y had to search every feasible 
integer point in y space. As k becomes larger, however, 
the number of enumerable points in y space becomes 
smaller for a given d and the problem can be solved 
more quickly. 

Doing the remaining substitution for w in PI 
gives the following form of the CPI-Block capacity 
model. 

P7: min Z = E [p^x.+h.v. ] + -[ E c'y +K] - iiSzilc,w 

7 1=1 1111 ic 1^1 11 k 1 o 



s . t . 

(P3.1), (P3.2), (3.1), and (3.2) 

X, V - 0 and INTEGER 
y > 0 

where c^ and K are as in P3. The basic algorithm is 
valid for P7 with a few changes. In step 4b let 
yj ■ yj "^he global upper bounds are changed 

slightly by the parameter k, since the maximum value 



for w^ is now 
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So using (3.2), (2.5) of Chapter II becomes 




N 



In order to satisfy (3«1)» (3«2) and y - 0, the lower 



bounds on y are first obtained from (3.1). If y - 0, 
then (3.2) Is also satisfied. If y^^ < 0, 1 - 1 - N, 
then by adding to y^^ for successive values of 
n^ = 1,2,3,..., until y^^ - 0, a feasible lower bound 
can be obtained which satisfies (3.2), 

CIDPI~Block capacity model 

Now for the CIDPI model, (2,11) can be re- 
written In terms of block capacity as 



1 

(w^-r^^) = -[d^ ^1 ” 

k 



1 



(w2-r2) « -[d2 - d^ - 2y^ + y2 + '^o^ 



(wi-ri) - i[d^ 



"^1-1 + ^1-2 - 2yi-i + yi^ 

1=3 



, , , , , 



N 



Substituting for r - 0 In P2 gives 




and INTEGER 
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Wjj — — 1 1 ^ ^N— 2^ ** ^ 3.nd INTEGER 

Ic 

Following a development similar to the first part of 
this section, the constraint is equivalent to the 

two constraints 

and 

Yj » njk + dj_i - dj + 2yj_i - yj_2 

where nj Is any Integer, but now nj < 0 Is allowable 
since the right-hand side of the first constraint Is 
no longer required to be a non-negative. 

In general, the two constraint sets which 

result are 

wi > i[d^+y;L-( w^+Vq)] 

^ 1 

(3.^) ^2 ~ i^l^^2“^i“2yi+y2+VQ] 

Jv 

and 

Yl = n^^k - dj^ + (kw^+Vo) 

(3.5) 

Y2 “ n2k + 2y^ + d^ - d2 - Vq 
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» n^k + 2 y 2 ~ yi ^2 ” ^3 INTEGER 

. 1 = 1,2,. . . ,N 

+ 2yN_i - yw_2 + dj^-i " 

Substitution for r in the remainder of P2 Just repeats 
the first three constraint sets of P5 and gives the 
model 

N 

P8: min Zq = ^[PiX^+hj^Vj^+f ”y^+( c^-f ^)Wj_ ] + M* 

s. t. 

(P5.1), (P5.2), (3.^), and (3-5) 

7, 7, w ^ 0 and INTEGER 
7 > 0 

where 

f = if* 
i k i 

M* = iw - liSzIlf w 

k k 1 o 

for f* and M from P5. 

1 

Modifications to the basic algorithm for P8 
are the same as for P7 except that, if the lower bound 
for any y^ is positive, then (n^k) should be subtracted 
from y^ for successive values n^^ » 1,2,3,..., until the 
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smallest - 0 is obtained. This corresponds to the 
unrestricted n in constraint set (3.5) which is caused 
by the possibility of having r^^ > 0 in any period. 

It should be noted here that the special 
cases of Chapter II are valid here only for k = 1. 

Another point worth noting is the increase in danger of 
rounding off from the related linear programming solu- 
tion as k increases. Not only are the unit capacity 
costs larger for large k, but also shifting a unit from 
one period to another has a proportionally greater 
effect on the optimal production-inventory schedule. 

Finally, the increase in the lattice Interval 
with k leads to a useful scheme when deciding on a capac- 
ity increase (and decrease) strategy. If several sizes 
of capacity are available, then the algorithm can be run 
for each value of k with the running time probably de- 
creasing for each larger value of k. This allows a third 
dimension of optimization for management with the cost 
of gaining the added Information governed by the algo- 
rithm running time for the smallest value of k. Note, 
however, that the form of Standard Transformation 1 does 
not permit k to be non-statlonary within a model. 

III. 3 Convex Production Cost 



Up until now, it has been assumed that all 



production was being carried out by one shift working 
regular hours with unit production costs p. In this 
section the CPI and CIDPI models with block capacity 
will be extended to allow multi-shift and overtime pro- 
duction. Although the work done here is for one addi- 
tional shift, further extensions are trivial, as will 
be pointed out later. 

The crucial assumption of this section is that 
the unit cost of production increases for each shift 
and for overtime. In the manufacturing context, this 
assumption is not restrictive and actually models real- 
ity, since it is commonly accepted that there exist 
"shift differentials" in manufacturing costs for a 
second and third shift. These can be ascribed, in part, 
to the increased time spent on maintenance and to 
greater human error in succeeding shifts. As for over- 
time production, the wage differential alone is suffi- 
cient to Justify the assumption. 

Let 

x|- second shift production in period 1 
p^- second shift unit production cost in 
period i 

Production capacity has been defined in terms of regu- 
lar time production. It is assumed that, in terms of 
the second shift, an equal amount of capacity is 
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available. Thus 




1 » 1 



Assume that 




i » 1 



» • • • 



Then the production cost function for period i 



is plecewise-linear convex. Similarly, with a suitable 
rearrangement of terms, total production cost 



is piecewise-linear convex. 



CPI^Block capacity- Convex production cost model 



First the extension of the CPI model with 

N 

block capacity will be presented. The term E p'x' 

1«1 i i 

must be added to the objective function of P^^, while 
the following constraints are added to the constraint 




N 



P = Z P.(x^,x*) 

iol 1 1 i 



set: 



t 
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i . 

X* - k Z w, - kw 1 = 

1 J»1 J ° 

(PI. la) 

X' - 0 and INTEGER 

and the appropriate block capacity coefficients are in- 
serted in the remainder of the model. Then (Pl.l), 

(PI. la), and (PI. 2) can be combined to give a set of 
necessary and sufficient conditions for feasibility 
which take into account two-shift production. 

2kw^ - Vi ^ dj^ - Vq - 2kw^ 

2k(2w^+W5) - V- - d, + d, - V. - Mkw 
1^2120 o 

• 

> N 

2k(Nw +. . .+w„) - V - Z d. - V - 2Nkw 
1 N N 1 o o 

Again by Standard Transformation 1 a solution for w can 
be obtained. Then substituting into w - 0 of PI gives 
the following two sets of constraints 

S'! - ''o + 

(3.6) yj i 2y^ + <^1 - <*2 - ■'o 

yj - 2yi_j^ - yi_2 + i - 3.....N 



and 
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= 2n^k - + 2 KWq + 



72 “ 2n2k + 2 ) 1 ^ " ^2 ” '^o 



(3.7) 



ri4 - 0 

and INTEGER 

i = 1, 2, . . . ,N 



= 2n3k + 2^2 ~ ■*■ ^2 ” ^3 



> • • • > 



^ ’ yN-2 ^ ^N-1 - 



Prom (3*7) it is clear that the lattice Interval on the 
y space has been increased to 2k with the Introduction 
of a second shift. 

Completing the substitution and transformation 
from PI gives 



P9 




s. t . 






1*2 



N 




i » 2 



N 



(P9.3) * ^1 




(3,6) 




I 
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(3.7) 

X, x', V - 0 and INTEGER 
y - 0 



where 



and 




i « 1, 



,N 



K'= ^ K - w 

^ ^ 2k ^ l'* o 



for c’ and K from P3. 

1 

The addition of the extra term to the ob- 
jective function and the extra set of constraints 
modifies the transportation problem structure of the 
subproblem. In effect, N additional sources have been 
added and the number of arcs has been doubled. This 
development might make it beneficial to reformulate 
the sub-problem as a transshipment problem which can 
be solved with an appropriate minimum cost flow algor- 
ithm (Ford-Fulkerson 62). An investigation would have 
to be made of the tradeoff between the faster algorithm 
but higher network growth rate of the tremsportation 
problem against the slower algorithm but lower network 

growth rate of the transshipment problem. Other algo- 
rithm changes and remarks on the bounds for y are 



identical to those for the CPI model of the previous 
section. 
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CIDPI - Block capacity - Convex production cost model 



The convex production cost case of the CIDPI 

model with block capacity can be obtained easily by a 

similar development. In P2, the terra ? p'x' is added 

i=l i i 

to the objective function, the constraints 



(P2,la) 



x' - k Z (w.-r.) - kw i = 1,...,N 

i j=l J J o 

x’ - 0 and INTEGER 



are added to the constraint set and the appropriate 
block capacity coefficients are inserted in the remain- 
der of the model. A set of necessary and sufficient 
conditions for feasibility is then obtained from (P2.1), 
(P2.1a), and (P2.2) as 

2k(wi-ri) - v^ ^ d^ - Vq - 2kw^ 

2k[2(wj,-r^) + (w2-r2) ] - V 2 - d^ + d2 - - ^kw^ 

2k[NCw^-r^) + ... + (Wj^-rj^)] - - J^d^ " ~ ^^^o 

Using Standard Transformation 1 to solve for r and 
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substituting for r - 0 in P2, the following two sets 
of constraints, analogous to (3.^) and (3.5), are ob- 
tained. 



( 3 . 8 ) 



w - -lc£d,+y.-v -2kw ] 

1 2 l lo O'* 

w, - -k[d.-dT-2y.+v ] 

2 2 2 1 1 0 

i'*’yi“^yi-i'*'yi— 2^ ^ * 3 ,...,n 



and 



( 3 . 9 ) 



'N 



2n^k - 


d, + V + 2kw 
1 o o 




2n2k + 


2yj . <J^- dj - 




2n^k + 


2yj - - dj 


n^ INTEGER 



1 = 1,2,. ..,N 

2nj^k + 2yM_i - yj^_2 ^N-1 ~ 



Again note that the lattice Interval on the y space 
has been Increased to 2k. 

Completing the substitution for r In P2, the 
model is transformed to the following: 




PIO: min Z 
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s. t. 



(PlO.l) 



^ ^ ^1 - ''o' 

^ ^ ^1 - ^l-l' 



i = 2,. . . ,N 



(P10.2) 



>'1 - - ''o' 

'‘i ■ ' ^i-i' 



1 = 2,. . . ,N 



(P10.3) 



+ ^; + '^i-l - ^1 = ^1 



1 = 1,...,N 



(PIO. 4) 



+ y-|L 



d± + y± - Yi.i 



- 0 



1 = 2,...,N 



(3.8) 



(3.9) 



X , X * , V , w - 0 and INTEGER 



y 



0 



where 




1 = 1,. . . ,N 



and 
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. 1 . „ . 

where and M are from P5. 

The remarks of the previous section relevant 
to the changes In the algorithm and In the bounds apply 
here also. As noted there, the subproblem can be form- 
ulated as either a transportation problem or a trans- 
shipment problem. 

As has been mentioned, the techniques of this 
section easily extend the models to encompass a third 
shift and overtime work. The capacity available for 
overtime can be assumed to be equal to that for one 
regular time shift. This presents a problem since 
overtime usually Includes both tlme-and-one-half and 
double-time wage periods. Some subjective unit pro- 
duction cost would have to be chosen to represent the 
entire period rather than differentiating between the 
two cost levels. 

Finally, as was already remarked, an add- 
ition of another shift Increased the lattice Interval 
on the y space. If the extension Is continued, then 
the first term on the right-hand sides of (3.6) and 
(3.8) becomes 3nk with a third shift, and 4nk with 
overtime. Increasing the interval by nk each time. 



Chapter IV 



The basic Idea of the enumeration algorithms 
presented Is similar to work done by Murty (Murty 68) 
and Gray (Gray 71) • A discussion of these related 
approaches opens this chapter. Also Included are a 
short discussion of the restart capability of the 
algorithms and a presentation of computational results. 
Finally, some conclusions regarding the preceding work 
and suggestions for possible extensions are given. 

rV.l The Decomposition Approach to Integer Programming 

Since the algorithms developed In Chapter II 
use the Idea of decomposition. It would be useful at 
this point to put them In perspective relative to some 
similar approaches In Integer programming. This Is not 
done to claim extensions to existing theory, but rather 
to show contrasts and similarities occurlng when a 
common approach Is used on different problem structures. 

One approach which uses the properties of 
duality theory to develop a gradient algorithm for 
mixed Integer programming was given by Benders (Hu 69). 

A drawback to the algorithm Is that efficient use of the 
gradient Information requires the solution of a series 
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of Integer programs where the number of constraints Is 
Increasing by one at each iteration. The algorithms 
of Chapter II only make partial use of the available 
gradient Information In the convexity properties which 
are developed. Possible further use of this Information 
was not Investigated In the course of this work. 

Other approaches more closely related In 
principle and In spirit to the algorithms of this 
thesis are the work of Murty (Murty 68) and Gray (Gray 
71). Both works deal with an exact solution to the 
fixed charge problem and, as pointed out by Gray, the 
two have mutually enhancing properties and could pro- 
bably be combined to give a powerful algorithm. 

Murty develops a decomposition approach which 
ranks the extreme points of the convex hull according 
to the value of the variable cost. An upper bound on 
variable cost Is given by the use of a minimum value for 
the fixed cost In a manner similar to that which gives 
an upper bound on capacity cost In the algorithms of 
this thesis. The variable cost upper bound can then 
be used to terminate the enumeration of extreme points. 
The algorithm's efficiency Is greatest when the fixed 
costs are small both In absolute and relative magnitude 
with respect to the variable costs (Gray 71). 

Up until now, no mention has been made of the 
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fixed charge aspects of PI and P2. It is not hard to 
make such an Interpretation of the capacity costs how- 
ever, since once Incurred, they are spread over the 
units produced In all future periods. The capacity 
cost of this thesis can thus be thought of In terms of 
a fixed charge although the problem formulation does 
not make this explicit. 

In light of the above, the approach of this 
thesis Is much closer to that of Gray than of Murty, 
since Gray's decomposition approach works mainly on 
tne fixed cost section of the problem and deals with 
the variable cost as a subproblem. His fixed cost 
section Is concerned with the optimal profile of 
open and closed routes of a transportation problem 
where each open route Incurs a fixed charge. The scheme 
used Is enumeration In a number system of base 2*^ (where 
n Is the number of sinks) which takes advantage of both 
the transportation problem structure and the 0-1 prop- 
erty of the Integer variables. In Gray's approach, an 
upper bound Is developed on the fixed cost which Is 
analogous to a of Chapter II. 

The major difference between Gray's (and 
Murty 's) model and those of this thesis Is the 0-1 
nature of the Integer variables. Transformation of PI 
and P2 to 0-1 variables could be done, but at the 
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expense of a huge Increase In the number of variables. 

Gray also encounters another difficulty not 
present In this thesis. He has no means of ensuring 
that each node of the enumeration Is feasible to the 
to the main problem, thus expending resources In enumer- 
ating, and testing for. Infeasible nodes. As has been 
mentioned previously, enumeration on y assures feasibil- 
ity at every node In the algorithms of this thesis. 

rv.2 Restart Capability of the Algorithm 

The algorithms developed have a feature that 
enhances their flexibility and usefulness, that Is the 
capability of starting them at any point In the enumer- 
ation. This feature Is Inherent in the method of 
enumerating on y; since given any y* and the value of 
the algorithm pointer (PNT) plus some value for the 
current bound on capacity cost (a; and the current op- 
timal cost, the algorithm can be started from Step 5 
and run normally. Even If a and current optimal cost 
are set to Infinity, the next execution of Steps 6 
through 9 will establish a new set of. bounds, though If 
the optimal solution precedes the point at which the 
algorithm Is restarted, the new bounds may not be as 
tight as those which would have resulted from the 
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optimal solution. If the optimal solution Is given by 
some vector after y It will be found whether or not 
Initial bounds are provided. 

This restart capability, when combined with 
the linear programming bound, provides a very useful 
tool to management, especially when a large problem 
which could be expected to have a long running time Is 
encountered. Instead of starting the algorithm from the 
lower bounds of y, a y* corresponding to a w slightly 
lower than any rounded-off solution to the linear pro- 
gram can be chosen as a starting point. The algorithm 
can then be run for a sufficient period of time to In- 
spect all w which could have been rounded-off from the 
linear programming optimum. An optimal solution Is 

, If none such Is 
found, then a management decision can be made based on 
the difference between the current optimum and Z* as 



assured here If It equals Z* 

' LP 



r- 



LP 

to whether the full algorithm should be run In an 

attempt to achieve further Improvement. 

It might also be advantageous, e.g. because 

of limited computer time, to run the algorithm for 

relatively short periods of time. If no terminal 

condition has been reached then, again, a management 

decision can be made on whether to attempt to Improve 

the current optimal solution with respect to Z* . 

LP 
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If so, the algorithm can be started right where It last 
terminated with no loss of efficiency. 

IV . 3 Employment Smoothing 

The relevance of the literature on employment 
smoothing was mentioned In Chapter I. Now that the de- 
velopment has been completed In the context of capacity 
expansion, a reinterpretation of the results can be made 
to handle the employment smoothing problem (O'Malley 
et.al. 66, Sliver 67, Llppman et.al 6?a and 67b, Llppman- 
Yuan 69 )-. Reference, In this section, will be to the 
most general models developed, l.e. those of section 
III. 3, from which all special cases can follow trivially. 

The time period for the employment smoothing 
problem Is much shorter than that for capacity expan- 
sion by the simple fact that hiring, firing, and usually 
even training, labor does not take as long as changing 
fixed capacity. It Is assumed, therefore, that fixed 
capital capacity Is never a constraint on the size of 
the work force. 

The work force units can be either In terms 
of Individuals, or groups such as gangs, crews, etc. 

It Is assumed that there Is a capacity associated with 
each labor unit which can be expressed In terms of the 



109 



consumptive factor bundles This capacity could be 

gotten from the associated fixed capital optimal oper- 
ating level capacity or by a separate marginal pro- 
ductivity level argument for labor as was done for 
capital In Chapter I. It Is again assumed that the 
production function exhibits constant returns to 
scale In labor and the consumptive factor. The accom- 
panying assumption Is that the existing labor force can 
be used at any level up to Its capacity with no lost 
opportunity cost of Idle labor other than wages. 

A change In the factors of capacity will have 
to be made. Before, the capacity definition was ex- 
pressed In terms of one shift doing regular time work. 
Now capacity Is defined In terms of all possible regular 
time output which can be achieved by the existing work 
force. Thus It Is assumed that the shifts are filled 
sequentially, l.e. before adding second shift labor, a 
full first shift must be operating. 

The CPI model P9 would have limited appli- 
cation In the employment smoothing sense because It 
Ignores capacity reduction. Its simplifying features 
relative to the CIDPI model however, make Its use at- 
tractive when the circumstances permit. Two Instances 
when such use would be feasible are those of a firm 
phasing In a new production line where the phase-ln 
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period Is relatively long, say a year; and that of a 
firm, with a reduced work force on an existing line, 
facing a projected growth In demand. 

Reinterpretation of P9 to employment smoothing 
leaves the Inventory variables v and their coefficients 
h unchanged. The production variables x now represent 
all regular time production, and no differentiation In 
production costs p among shifts Is allowed. Overtime 
production Is x’ and again undifferentiated overtime 
production costs p’ must be assigned. The capacity ex- 
pansion variables w now become work force hiring var- 
iables and the associated costs c can Include all the 
normal costs associated with hiring personnel, l.e. ad- 
vertising, personnel staff costs, training, etc. Wages 
now take the place of maintenance costs for m. 

As the reader has probably noticed, an Inter- 
pretatlonal difficulty exists concerning the capacity 
ascribed to overtime production In P9. The assumption 
for that model was that overtime capacity was approx- 
imately equal to the capacity of a regular time shift. 
What Is necessary now. If Standard Transformation 1 Is 
to be used. Is the seemingly ludicrous assumption that 
overtime capacity Is equal to all regular time capacity. 
With a time span of one to one and one-half years, how- 
ever, discounting will not be significant so that over- 
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time production costs will be much higher than any reg- 
ular time costs over the period. Furthermore, hiring 
costs are now In the same order of magnitude as pro- 
duction costs so that, within a reasonable range, 
hiring additional regular time labor should be cheaper 
than working any significant amount of overtime. A 
reasonable approach to the problem would be to first 
solve It for regular time only, and then try a solution 
permitting overtime. If the overtime assumption causes 
difficulty, then either a subjective modification of 
costs can be made and the problem re-solved or the 
approach has to be abandoned. For the majority of 
cases It would seem that the approach Is useful. 

For the more realistic CIDPI model PIO, the 
only additional Interpretation necessary Is that the r 
now represent the firing variables and the f the firing 
costs. The same difficulties exist for overtime capacl 
ty and the previous remarks apply equivalently here. 

An argument can be made for the equivalence 
of hiring and firing costs In any period. Certainly 
what Is Included a..d excluded from each cost Is sub- 
jective, Assuming that equality exists and, further, 
that the c are monotone non- Increasing with time, then 
the special case of the network formulation Is appli- 
cable to the employment smoothing problem with the 
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preceding interpretations on the relevant variables 
where in addition w^^ < 0 corresponds to the firing 
variable r^^^ > 0 of PIO. Moreover, use of the network 
excludes the possibility of overtime and requires that 
k ® 1. No obvious extension of the network model to 
remove these restrictions exists. 

Finally, the assignment of a value for k in 
P9 and PIO will undoubtedly be more difficult in terms 
of labor than it was in terms of fixed capital. On the 
other hand, the sensitivity of the models to k should 
not be significant for labor because the hiring and 
firing costs are much smaller than the corresponding 
capacity change costs, and thus the Impact of fluctu- 
ations in the labor force level will not be as great. 

This work, then, becomes the first exact 
Integer approach for the employment smoothing problem.^ 
The algorithms given for the special demand cases in 
sections II. 6 and II. 9 are more general than those of 
Lippman et.al. (Llppman et.al. 67b), and algorithms to 
handle general demands and general linear costs have 



^As was mentioned in Chapter I, a "round-off" 
approach from a linear programming model has been pre- 
sented (O’Malley et. al. 66). 
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been given in addition. Finally, the similarity of 
the two areas of capacity expansion and employment 
smoothing certainly suggests that a combined effort 
would be useful. 

IV. 4 Computational Results 

The algorithms developed In this thesis were 
coded In Fortran IV and run on the IBM 7094 In the 
Johns Hopkins University computing center. The test 
problems used were set up by the author using non-random 
data. Twelve problems were generated In all, each one 
as a ten-period problem. When run, however, each prob- 
lem was run eight times using the first k periods of 
data each time for k = 3,..., 10. Thus e.g. for problem 
10, the three period problem Is numbered 10.1, while the 
ten period problem Is numbered 10.8. This numbering 
system Is followed in the tables of data which follow. 

The following conditions are satisfied by all 
of the problems: a) Costs are monotone non-lncreaslng 
with time. In addition, the costs of Problems 7* 9 

and 10 are stationary, b) k = 1 (1 Capacity unit = 1 
demand unit), c) Only one shift of regular time work Is 
being used. 

A linear program was solved for the CPI and 
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CIDPI case of each problem (using models PI and P2). 
Because of the •'nearly” unlmodular structure of the 
general models, a large number of Integer extreme points 
probably exist In each convex hull. In fact, for the 
166 linear programs solved for Problems 1 and 3 to 12, 
a total of 40 terminated Integer. 

It had been thought that the number of periods 
would be the most significant factor In problem running 
time and that the appropriate curves could be presented. 
Instead, the models turn out to be sensitive to both the 
cost and demand structures of the Individual problems, 
as well as being dependent on the number of periods. 

The computational results are, therefore, presented In 
tabular form In Tables IV-1, IV-2, and IV-3 at the end 
of this section. 

An example of the sensitivity of the algorithms 
to demand and cost structure Is given by Problems 9 and 
12 using the CPI algorithm. The data for the problems 
Is as follows: 

Problem £ 

d = (5, 7, 11, 6, 10, lA, 9, 13, 19, 16) 
c = ( 30, 30, . . ., 30) 

ffi = ( 3, 3, ..., 3) 

P = ( 5, 5, . . . , 5) 

E = ( 1, 1, ..., 1) 
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Problem 12 

5 - (2, 3. 5, 4, 6, 8, 7, 9. 12. 10) 
c - (250, 240, 230, 220, 210, 200, 190, 

180 , 170 , 160 ) 

m = (3, 3, ...» 3) 

P = (5, 5, 5, 4, 4, 4, 3. 3, 3) 

h = (1, 1, .... 1) 

For t = running time 

n = number of periods 

a curve fitting program was run for a polynomial 

t = an*^ 

and an exponential 

. bn 

t = ae 

fit on the data. For Problem 9» using 7 data points the 
results were: 

Polynomial 

a = .00028 b = 5.10 
Coefficient of correlation = .960 

Exponential 

a = .00635 b = .94 

Coefficient of correlation = .606 
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while using 7 data points for Problem 12 the results 
were: 



Polynomial 

a * .00006 b » 5.63 
Coefficient of correlation = .875 

Exponential 

a * .00573 b = .87 
Coefficient of correlation = .973 



When the data for 6 common points was averaged the 
following resulted: 

Polynomial 

a = .00003 b = 6.21 

Coefficient of correlation = . 561 
Exponential 

a = .00339 b = .993 
Coefficient of correlation = .770 

Similar results were obtained using the data from the 
other problems for both the CPI and CIDPI cases. 

The only conclusion which can be drawn. 
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then. Is that the number of periods involved is a sig- 
nificant factor in problem running time, but there 
exist other significant factors. Inspection of the 
tables will verify the claim made at the end of section 
II. 8 that longer running times should result for the 
CIDPI algorithm. 

Two interesting examples remain to be shown, 
the first demonstrating the acknowledged fallacy of 
rounding-off the linear programming solution, and the 
second a demonstration of the power. of the restart ca- 
pability discussed in section IV. 2. 

Example 1: For the CPI case of Problem 12.2, the vari- 

able values resulting from the linear and Integer pro- 
gramming optima respectively are (with w^ = 1 and Vq = 1) 

LP ® 2.25* W2 = w^ = w^^ = 0 

= X2 = = x^ = 3.25 

v^ = 2.25, v^ = 2.5, = .75, Vi|= 0 

0, = 3, w^ = w^ = 0 

1, X2 = X3 = X4 » 4 

0, ''^2 ~ 1> V3 = vij = 0 

The difference between building all capacity in period 
1 and building it all in period 2 could be significant 



IP w « 
— 1 

^1 = 

vi = 
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to management. Furthermore, the product ion- inventory 
policies are quite different in each case. 

Example The CPI algorithm for Problem 8.8 took 
27.067 seconds to terminate where termination resulted 
from the linear programming bound. The optimal capacity 
expansion vector for the linear programming solution 
was 

w = ( 1 , 4 , 3.875, 0, ..., 0 ) 

The CPI algorithm was started with a y corresponding to 
w-j^ = 1, W 2 = 4, w^ = 3, and the remaining y^ through y^^Q 
set at their lower bounds. The time to termination by 
the linear programming bound from this starting point 
was 5.9167 seconds. 

In terms of core usage in the computer, the 
algorithms are limited only by the requirements of the 
transportation problem. The version of the transporta- 
tion algorithm used requires only one mxn array where 
for the algorithms of Chapter II, m = n = N, while if 
a transportation algorithm is used for the convex pro- 
duction cost use of section III. 3 , m ** 2 N, and n ® N. 

The enumeration itself is performed using a single 
vector N +2 words in length. The growth in core require- 
ments is therefore dominated entirely by the growth of 
the transportation array. 
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Table IV-1 



Computational results for the CPI algorithm 
with monotone non-decreasing demand 



Problem 


Periods 


Optimal Values 


Run Tlme(sec) 


Step 8 


lumber 


(N) 


LP 


IP 


LP 


IP 


’Termination 


1.2 




385.0 


388 


.73 


.383 


No 


6.2 


H 


176.0 


179 


.77 


.20 


No 


7.2 


H 


211.0 


217 


.77 


.20 


No 


1.3 


5 


557.0 


558 


1.05 


1.533 


No 


6.3 


5 


250.0 


252 


1.03 


.40 


No 


7.3 


5 


309.0 


317 


1.05 


.317 


No 


l.ii 


6 


735.25 


736 


1.33 


.733 


Yes 




6 


335.25 


336 


1.33 


.283 


Yes 


7.^ 


6 


429.50 


437 


1.47 


.717 


No 


1.5 


7 


884.60 


885 


1.83 


1.383 


Yes 


6.5 


7 


418.40 


420 


1.82 


3.083 


No 


7.5 ' 


7 


572.85 


578 


2.35 


2.4 67 


No 


1.6 


8 


1052.0 


1055» 


2.27 


8'*’ min 




7.6 


8 


739.75 


741 


3.02 


11.383 


No 


6.7 


9 


583.60 


586 


3.50 


85.233 


No 


7.7 


9 


940.75 


942 


3.75 


75.350 


No 


6.8 


10 


668.25 


670» 


4.23 


ll'^ min 




7.8 


10 


1173.0 


1174 


4.67 


474.50 


No 



*Thls was the current optimal solution when the 



enumeration was stopped. 

^Step 8 uses r^^pi 



as a termination rule for the 



enumeration. 
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Table IV- 2 

Computational results for the CPI algorithm 
with general demand 



Problem 


Periods 


Optimal Values 


Run Tlme(sec) 


Step 8 


Number ' 


(N) 


LP 


IP 


LP 


IP 


Termination 


4.1 


3 


299.33 


310 


.40 


.283 


No 


9.1 


3 


373.0 


391 


.40 


.20 


No 


4.2 


4 


472.0 


478 


.73 


.883 


No 


8.2 


4 


271.0 


277 


.75 


.317 


No 


9.2 


4 


425.0 


445 


.72 


.383 


No 


12.2 


4 


668.75 


821 


.75 


.250 


No 


4.3 


5 


578.0 


581 


1.07 


2.533 


No 


8.3 


5 


377.0 


385 


1.07 


.70 


No 


9.3 


5 


513.0 


521 


1.10 


.517 


No 


10.3 


5 


232.0 


240 


1.07 


.433 


No 


12.3 


5 


855.40 


893 


1.07 


.450 


No 


4.4 


6 


636.0 


636 


1.43 


2.283 


Yes 


9.4 


6 


673.0 


678 


1.50 


1.650 


No 


10.4 


6 


313.50 


327 


1.52 


1.0 


No 


12.4 


6 


1099.0 


1185 


1.52 


1.050 


No 


8.5 


7 


568.14 


572 


2.35 


3.517 


No 


9.5 


7 


747.43 


750 


2.40 


2.10 


No 


10.5 


7 


346.0 


362 


2.30 


1.267 


No 


12.5 


7 


1245.86 


1268 


2.37 


1.850 


No 



I 



j 
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Table IV-2 (cont.) 



Number 


Periods 


LP 


IP 


LP 


IP 


Step 8 


4.6 


8 


781.0 


786 


2.87 


35.583 


No 


8.6 


8 


734.0 


735 


3.05 


19.60 


No 


9.6 


8 


882.50 


884 


3.12 


3.867 


No 


10.6 


8 


421.0 


422 


3.08 


1.817 


No 


12.6 


8 


1449.0 


1534 


3.05 


4.233 


No 


8.7 


9 


843.25 


844 


3.83 


4.550 


Yes 


9.7 


9 


1101.25 


1102 


3.83 


16.433 


Yes 


10.7 


9 


523.0 


525 


3.80 


7.20 


No 


12.7 


9 


1738.89 


1844 


3.90 


21.867 


No 


4.8 


10 


921.75 


929* 


5.43 


10'*' min 


No 


8.8 


10 


1022.75 


1023 


5.62 


27.067 


Yes 


9.8 


10 


1261.50 


1263 


5.48 


287.550 


No 


10.8 


10 


603.25 


604 


5.45 


3.817 


Yes 


12.8 


10 


1917.50 


1965 


5.67 


41.50 


No 
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Table IV- 3 

Computational results for the CIDPI algorithm 
with general demand 



Problem 


Periods 


Optimal Values 


’■R”n"T3 


.me ( sec ) 


Step 8 


Number 


(N) 


LP 


IP 


LP 


IP 


Termination 


4.1 


3 


299.33 


302 


.57 


1.067 


No 


9.1 


3 


373.0 


387 


.57 


.633 


No 


C\J 

• 


4 


472.0 


472 


1.12 


10.80 


Yes 


8.2 


4 


271.0 


276 


1.13 


4.650 


No 


9.2 


4 


415.67 


429 


1.13 


9.217 


No 


12.2 


4 


668.75 


728 


1.15 


2.267 


No 


4.3 


5 


578.0 


581* 


1.73 


9'*’ min 




9.3 


5 


513.0 


521 


1.77 


255.350 


No 


10.3 


5 


232.0 


240 


1.73 


28.167 


No 


12.3 


5 


855.4 


868 


1.73 


38.350 


No 


9.4 


6 


673.0 


69 5» 


2.55 


8‘*’ min 




10.4 


6 


313.5 


327 


2.53 


616.980 


No 


12.4 


6 


1099.0 


1113 


2.53 


880.0 


No 
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IV . 5 Conclusions 



The algorithms developed in this thesis take 
advantage of the structure inherent in the combined 
capacity expansion and production-inventory problems 
when no backordering Is allowed. The decomposition 
approach of handling the fixed-charge (capacity ex- 
pansion) and variable cost (production-inventory) por- 
tions of the problem separately allowed the application 
of special techniques for each part which could then 
be combined into an efficient algorithm. The regular 
behavior of the y variables in the fixed cost por- 
tion of the models allowed the development of useful 
optimal cost convexity properties in the variable cost 
section and in the overall problem. 

The main problem of Integer programming al- 
gorithms has not, unfortunately, been avoided in this 
work. That is, these algorithms display the same 
unpredictability in running time displayed by all 
known Integer algorithms. Furthermore, computational 
experience has shown that the algorithm running times 
are sensitive to demand and cost structures as well as 
being dependent on the number of periods involved. 

The flexibility provided by the restart 
capability of the algorithms is one of their most 
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attractive features from a management point of view. 

The feature allows for efficient use of the linear 
programming solution initially, by permitting a 
search among all rounded-off values. Available com- 
puter time can then be used efficiently for the main 
algorithm by running the program for short periods of 
time, since no information is lost when the enumeration 
is stopped. 

The economic problem of lumpiness in capacity 
acquisition has been handled rather easily in the 
extensions of Chapter III. In fact, the more severe 
the lumpiness is, l.e. the higher the value of k, the 
larger the Intervals in y-space where the enumeration 
is performed. 

Finally, the observation of an equivalence 
between the CPI and CIDPI models and the employment 
smoothing literature is of no small significance. The 
amount of work done in employment smoothing has been 
fairly extensive and unification of the results of 
these two areas would enhance the existing results in 
both . 

IV . 6 Ext enslons 



The restrictive assumptions which led to the 
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results of this thesis provide for some obvious ex- 
tensions. Perhaps the one with the greatest usefulness 
in the management context would be a provision for 
backordering, where the backordering horizon could be 
restricted to M periods, for Me{ 1 , 2 , 3 , . . . } . Unfortu- 
nately, Incorporation of backordering in the model will 
destroy the structure of Standard Transformation 1 and, 
in addition, if a decomposition approach is taken to 
solve the problem, the size of the enumeration space 
would be greatly enlarged. On the surface it seems as 
if an entirely new approach must be developed for this 
extension. 

Another useful extension would be the develop- 
ment of a model which could handle concave costs of 
capacity expansion, l.e. decreasing marginal cost in 
the addition of capacity in any period. The optimal 
policy results of Manne and Velnott for capacity ex- 
pansion with concave costs (Manne-Veino tt 67) might 
be useful here. 

As was mentioned in Chapter I, the work in 
capital budgeting is Just one step up in aggregation 
from the content of this thesis. A simple extension 
in the context of the work already done would be to 
have a capacity expansion budget in each period. This 
adds N budget constraints to the models and these 
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constraints could be handled easily in the fixed charge 
(capacity) part of the enumeration algorithms. 

An extension to many production lines with an 
overall capacity budget in each period would not 
be so easy, even if each line had a single product. 

The difficulty will spring from the problem of opti- 
mally allocating the capacity resources among lines in 
a period as well as optimizing for each line over N 
periods. A very large-scale integer program will cer- 
tainly result. 

A smaller scale, but equally interesting, 
problem would be an extension to multiple products on 
the same production line. As was noted in Chapter III, 
the work of this thesis limits k to a single value for 
any model because of the form of Standard Transformation 
1. Multiple products are sure to imply different pro- 
duction rates and, hence, different values of k. The 
challenge appears interesting. 



127 . 



BIBLIOGRAPHY 

Bowman, E,H. (1956), "Production Scheduling by the 
Transportation Method of Linear Programming," Operations 
Research , vol.M, no.l, pp, 100-103 

Carroll, Stephen J, (I 968 ), "Optimal Production, Inven- 
tory, and Capital Accumulation Policies for the Firm," 
Unpublished Ph,D. Dissertation, The Johns Hopkins 
University, Baltimore. 

Cobian, M.J, and P.T, Lele (1970), "Long Range Planning 
of Multi-Plant Capacity," Proceedings of the 21st AIIE 
Conference , pp, 97-101 

Fetter, Robert B, (I 96 I), "A Linear Programming Model 
for Long-Range Capacity Planning," Management Science , 
vol. 7, no. pp. 372-378 

Ford, L.R, Jr, and D.R. Fulkerson (I 962 ), Flows in 
Networks , Princeton University Press, Princeton, N.J, 

Gray, Paul (1971), "Exact Solution of the Fixed-Charge 
Transportation Problem," Operations Research , vol. 19, 
no. 6, pp. 1529-1538 



128 . 



Hinomoto, Hirohlde (1965), "Capacity Expansion with 
Facilities Under Technological Improvement," Management 
Science , vol. 11, no. 5, pp. 581-592 

Hlrshlelfer, J. (1970), Investment . Interest . and 
Capital '. Prentice-Hall, Englewood Cliffs, N.J. 

Howard, G.T. and G.L. Nemhauser (1968), "Optimal 
Capacity Expansion,” Naval Research Logistics Quarterly . 
vol. 15, no. 4, pp. 535-550 

Hu, T.C. (1969), Integer Programming and Network Flows . 
Addlson-Wesley Publishing Company, Reading, Mass. 

Iglehart, Donald (1965), "Capital Accumulation and 
Production for the Firm: Optimal Policies,” Management 
Science, vol. 12, no. 3, pp. 193-205 

Kalro, A.H. and S.R. Arora (1970) , "Optimal Capacity 
Decisions for Production Facilities”, AIIE Transactions , 
vol. 2, no. 4, p. 321 

Llppman, Steven A., Alan J. Rolfe, Harvey M. Wagner, and 
John S.C. Yuan (1967a), "Optimal Production Scheduling 
and Employment Smoothing with Deterministic Demands," 



129 . 



MzinaRement Science , vol, 1^, no. 3, PP. 127-158 

Llppman, Steven A., , , and (1967b), 

“Algorithms for Optimal Production Scheduling and 
Employment Smoothing," Operations Research . vol, 15, 
no. 6, pp. 1011-1029 

Llppman, Steven A. and John S.C, Yuan (I969), 

"Discounted Production Scheduling and Employment 
Smoothing," Naval Research Logistics Quarterly , vol. 16, 
no. 1, pp. 93-110 

Manne , Alan S. (196I), "Capacity Expansion and 
Probablistlc Growth," Econometrlca . vol. 29, no. 
pp. 632-643 

Manne, Alan S. and Arthur F. Veinott, Jr. (1967), 
"Optimal Plant Size with Arbitrary Increasing Time 
Paths of Demand," Chapter 11 in Manne, Alan S. , 
Investments for Capacity Expansion ; Size . Location , and 
Time-Phasing , The M.I.T. Press, Cambridge, Mass., 1967 

Maxwell, W. David (I 965 ), "Short-Run Returns to Scale 
and the Production of Services," Southern Economic 
Journal , vol. 32, no. 1, pp. l-l4 



130 . 



Montogomery, Douglas C. (1971), "Stochastic Capacity 
Decision Models for Production Facilities," Proceedings 
of the 22nd AIIE Conference , pp. 315-319 

Murty, Katta G. (1968), "Solving the Fixed Charge 
Problem by Ranking the Extreme Points," Operations 
Research . vol, l6, no. 2, pp. 268-279 

O'Malley, Richard L. , Salah E. Elmaghraby, and John W. 
Jeske, Jr. (I 966 ), "An Operational System for Smoothing 
Batch-Type Production," Management Science , vol. 12, 
no. 10 , pp. B-il43-B-i|il9 

Sliver, Edward A. (1967), ”A Tutorial on Production 
Smoothing and Work Force Balancing," Operations 
Research , vol. 15, no. 6, pp. 985-1010 

Smith, Vernon L. (I 96 I), Investment eind Production , 
Harvard University Press, Cambridge, Mass. 

Velnott, Arthur F. Jr, and Harvey M. Wagner (1962), 
"Optimal Capacity Scheduling- Parts I and II," 
Operations Research . vol. 10, no. pp. 518-5^6 



Welngartner, H. Martin ( 1963 ), Mathematical Programming 



131 . 

and the Analysis of Capital Budgeting Problems , 
Prentice-Hall, Englewood Cliffs, N.J. 

Zabel, Edward (I963), "Efficient Accumulation of 
Capital for the Firm," Econometrlca , vol. 31 » no. 1 - 2 , 
pp. 131-150 



I 



132 . 



VITA 

Richard H. Barchl, a resident of Yardley, 
Pennsylvania, was born on 8 November 19^2 in Washington, 
D.C. His undergraduate work was done at the United 
States Naval Academy where he was awarded the B.S. in 
1966, His final year at the Naval Academy was spent in 
research in Fluid Mechanics on the topic of Boundary 
Layer Turbulence and Drag Reducing Additives.^ 

After graduation, he served two years as 
Missile Officer in USS Joseph Strauss (DDG-I6) stationed 
at Pearl Harbor, Hawaii. Matriculating to the Depart- 
ment of Operations Research of The Johns Hopkins 
University in 1968, the author completed the require- 
ments for the Ph. D. in March of 1972, and returned to 
duty at sea with the U.S. Navy. 



^Johnson, Bruce and Richard H. Barchl "Effect of 
Drag-Reducing Additives on Boundary Layer Turbulence," 
Hydronautlcs » vol. 2, no. 3, July, 1968, pp. 168-175. 
Johnson, Bruce and R.H. Barchl "Effects of Polymers on 
Boundary Layer Turbulence," Proceedin gs of the 11th 



International Towing Tank Conference , Tokyo, Japan 1967. 




I 



1 

I 



I 



I 



I 



I 



I 

I 



I 




Thes i s 
B2173 



132921 

Barchi 

An inteqer pro<^ran- 
ming approach to capa- 
city expansion and 
production planning. 



